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CONVERSION  FACTORS,  METRIC  (SI)  TO  U.  S.  CUSTOMARY  AND 
U.  S.  CUSTOMARY  TO  METRIC  (Si)  UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in  this  report  can  be  converted  as  follows: 


Multiply To  Obtain 

Metric  (SI)  to  U.  S.  Customary 


millimetres 

centimetres 

metres 

kilometres 

square  centimetres 

cubic  metres 

grams 

kilograms 

metric  tons 

grams  per  cubic  metre 

dynes 

centimetres  per 
second 


0.0394 
0.3937 
3.2808 
0.6213 
6. U516 
35.3147 
0.002205 
2.2051 
2205.0 
62.43 

2.248  x i(f6 
0.3937 


inches 

inches 

feet 

miles  (U.  S.  statute) 
square  inches 
cubic  feet 
pounds  (mass) 
pounds  (mass) 
pounds 

pounds  per  cubic  foot 
pounds  (force) 
inches  per  second 


metres  per  second 


3.2808 


feet  per  second 


U.  S.  Customary  to  Metric  (Si) 


feet  0.3048 

grains  0.0648 

pounds  (mass)  0.4535924 

pounds  per  square  inch  6.894757 

degrees  (angle)  0.01745329 

Fahrenheit  degrees  5/9 


metres 

grams 

kilograms 

kilopascals 

radians 

Celsius  degrees  or  Kelvins* 


* To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  read- 
ings, use  the  following  formula:  C = (5/9)(F  - 32).  To  obtain  Kelvin 

(K)  readings  use:  K = (5/9)<F  - 32)  + 273.15. 
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EXPERIMENTAL  VERIFICATION  OF  A THEORETICAL  LOADING 


FUNCTION  DESCRIBING  MOMENTUM  TRANSFER 
FROM  AN  EXPLOSION  TO  A TREE  STEM 


PART  I:  INTRODUCTION 

Background 

1.  An  important  characteristic  of  the  U.  S.  Army’s  tactical 
mobility  is  the  widespread  and  effective  use  of  troop-carrying  helicop- 
ters. Experience  gained  in  Southeast  Asia  indicated  that  the  Army  must 
develop  the  capability  to  conduct  effective  military  operations  in 
forested  areas  that  are  intensively  used  by  the  enemy  for  sanctuaries 
and  approach  and  retreat  routes.  Helicopter  landing  zones  (HLZ's)  are 
essential  in  such  operations  for  tactical  deployment  of  troops,  resup- 
ply, evacuation  of  wounded,  etc.  An  enemy’s  ingenious  use  of  various 
antihelicopter  mines  and  booby  traps  makes  the  use  of  existing  clearings 
for  HLZ's  dangerous;  furthermore,  establishing  an  HLZ  in  hostile  terri- 
tory by  conventional  engineering  techniques  is  both  hazardous  and  time- 
consuming.  Clearly,  a new  technique  is  needed  for  constructing  HLZ's 

in  forested  areas  that  would  meet  the  Army's  mission  requirements  with 
a minimum  exposure  of  men  and  equipment  to  a hostile  environment. 

2.  In  1968  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  began  a study  for  the  U.  S.  Air  Force  on  the  feasibility  of  using 
large  bombs  to  rapidly  construct  clearings  in  forested  areas.  Full- 
scale  experiments  conducted  at  Fort  Benning,  Georgia,  and  in  South 
Vietnam  demonstrated  the  utility  of  the  Air  Force  concept.  The  results 
of  these  experiments  were  used  to  structure  a mathematical  model  to 
predict  tree-remnant  height  as  a function  of  distance  from  the  center 
of  the  explosion,  weapon  yield,  tree  stem  diameter,  and  tree  strength. 
The  results  obtained  from  the  model  depend  largely  on  a theoretical 
loading  function  used  in  the  solution  of  a partial  differential  equa- 
tion that  describes  the  motion  of  an  individual  tree  stem  being  acteu 
upon  by  a shock  front.  Although  the  loading  function  was  structured 
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using  all  available  reliable  data,  the  initial  project  plan  did  not 
include  provisions  to  verify  that  the  loading  function  reliably  charac- 
terized the  action  of  the  shock  front  on  the  tree  stem;  however,  re- 
sults obtained  from  the  model  showed  good  agreement  with  test  results. 
The  output  of  the  model  can  be  used  to  compare  predicted  tree-remnant 
heights  with  vertical  ground  clearances  required  for  helicopter  landing 
and  takeoff.  This  comparison  will  indicate  whether  a clearing  is  satis- 
factory for  landing  zone  operations  prior  to  detonation  of  the  munition. 
The  development  of  the  model  and  applications  are  described  in  Refer- 
ence 1 and  Appendices  A and  B of  this  report. 

Purpose 

3.  The  purposes  of  the  study  reported  herein  were  to: 

a.  Develop  reliable  methods  for  evaluating  impulse  loading 
of  tree  stems  being  acted  upon  by  a shock  front. 

b^.  Acquire  adequate  data  to  validate  the  theoretical  loading 
function  used  in  the  partial  differential  equation  de- 
scribing tree  stem  motion  under  the  influence  of  a shock 
front or  if  the  loading  function  could  not  be  validated, 
use  the  data  to  generate  a new  loading  function. 

Scope 

1*.  Test  equipment  was  assembled  to  measure  the  impulse  loading 
of  a tree  stem  by  a shock  front.  Slash-pine  logs  were  used  to  repre- 
sent tree  stems.  Seven  explosives  of  0.09-metric-ton*  (200- lb)  TNT 

equivalent  yield  were  detonated.  Blast  data  from  each  shot  were  com- 

2-5 

pared  with  established  overpressure  and  time-of-arrival  (TOA)  data 
to  verify  that  the  desired  yield  had  been  obtained  and  that  the  shock 
front  propagated  away  from  ground  zero  (GZ)  in  a symmetrical  manner.  A 
method  was  developed  for  evaluating  the  theoretical  loading  i unction 
in  terms  of  the  impulse  experienced  by  the  logs.  Experimental  and 
theoretical  velocities  were  compared. 

* A table  of  factors  converting  metric  (SI)  units  of  measurement  to 
U.  S.  customary  units,  and  U.  S.  customary  to  metric  (Si)  is  given 
on  page  5. 
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PART  II:  EXPERIMENTAL  DESIGN  AND  PROCEDURES 


Test  Site 


5.  The  series  of  test  shots* **  was  conducted  on  a flat  plain  at 

the  WES  Big  Black  Test  Facility,  16  km  east-southeast  of  Vicksburg, 

Mississippi,  on  the  west  bank  of  the  Big  Black  River.  The  first  soil 

layer  at  the  test  site  generally  consists  of  a thin  mantle  of  topsoil 

underlain  by  silty  clay  (CL  according  to  the  Unified  Soil  Classifica- 

7 3 

tion  System)  with  a dry  density  of  2.3  x 10  g/m  . A seismic  refrac- 
tion survey  conducted  at  the  site  indicated  that  the  compression  wave 
velocity  in  the  first  layer  was  300  m/sec;  the  minimum  depth  of  the 
first  layer  was  6 m. 

6.  The  site  layout  for  a representative  test  shot  is  shown  in 
Figure  1.  The  center  of  the  explosive  charge,  hereafter  referred  to  as 
GZ,  was  75  m from  the  nearest  vertical  discontinuity  that  would  cause 
any  reflections  of  the  shock  front.  At  this  distance  the  overpressure 
had  decayed  to  less  than  1 psi.*# 

Explosive 

7.  Ammonium  nitrate  (mixed  with  6 percent  by  weight  No.  2 diesel 
fuel)  was  chosen  for  the  explosive  material  instead  of  TNT.  An  ammonium 
nitrate-fuel  oil  compound  (hereafter  referred  to  as  AN/FO)  is  oxygen- 
balanced  and,  hence,  the  chemical  reaction  goes  on  to  completion  at  the 
time  of  detonation.  A TNT  explosion  is  not  an  oxygen-balanced  reaction 
and,  therefore,  may  require  several  milliseconds  (msec)  to  attain  com- 
plete detonation  of  the  explosive  material.  Therefore,  use  of  AN/FO  im- 
proved the  probability  of  attaining  a very  rapid  rate  of  energy  release. 


* The  seven  test  shots  of  0.09-metric-ton  TNT  equivalent  yield  will  be 
referred  to  as  BB001-BB007  throughout  the  remainder  of  this  report 
when  it  is  necessary  to  specify  information  about  a particular  shot. 

**  psi  is  retained  as  the  unit  of  pressure  as  opposed  to  its  metric 
equivalent  because  this  unit  remains  in  common  use. 
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Figure  1.  Representative  layout  for  experimental  equipment 


Also,  the  safety  precautions  applicable  to  AN/FO  explosives  are  consid- 
erably less  restrictive  than  those  applicable  to  TNT.  The  AN/FO  explo- 
sive was  prepared  by  mixing  No.  2 diesel  fuel  oil  with  the  ammonium  ni- 
trate (Gulf  N-IV  fertilizer)  at  a ratio  of  6 parts  fuel  oil  by  weight  to 
94  parts  ammonium  nitrate.  The  ammonium  nitrate  was  sensitized  (oiled) 
at  the  site  by  vigorously  pouring  the  mix  back  and  forth  between  large 
metal  cans.  After  the  mixing  procedure,  the  explosive  was  stored  for 
at  least  24  hr  to  allow  maximum  oil  absorption  by  the  ammonium  nitrate.* 
The  average  equivalent  weight  of  TNT  to  AN/FO  in  the  40-psi  region  is 
O.85.  To  produce  a 0.09-metric-ton  equivalent  TNT  yield,  104.5  kg  of 
AN/FO  were  required. 

8.  For  test  shots  BB001-BB003  (see  Table  l)  the  AN/FO  was  placed 
in  a cubical  Styrofoam  container  (Figure  2a)  having  a spherical  cavity 
0.64  m in  diameter.  The  height  of  burse  (measured  from  the  ground  to 
the  center  of  the  cavity)  was  0.72  m.  This  value  was  scaled  down  from 
the  height  of  burst  of  the  Fort  Benning  shot,'*'  by  the  cube-root  scaling 
laws  applicable  to  explosives, ^ so  that  the  results  from  the  Fort  Ben- 
ning test  and  those  from  tests  conducted  at  the  Big  Black  Test  Facility 
could  be  compared.  The  explosive  was  detonated  by  a 4-lb  C-4  booster 
with  a No.  10  Engineer's  Special  Blasting  Cap  and  200-grain  Primacord. 
The  cubical  exterior  shape  of  the  Styrofoam  container  was  modified  to  a 
spherical  shape  (Figure  2b)  after  shot  BB003  (Table  l)  to  preclude  the 
possibility  that  the  additional  volume  of  the  cubical  container  might 
contribute  to  the  heat  of  the  explosion  or  could  in  some  way  distort 
the  symmetry  of  the  shock  front.  No  evidence  was  noted  on  subsequent 
tests  of  any  differences  in  the  explosions  as  a result  of  changing 
containers. 


* Gulf  N-IV  ammonium  nitrate  fertilizer  is  in  the  form  of  spheres 
(average  diameter  1.5  mm)  that  are  generally  called  "prills."  The 
spherical  surface  configuration  allows  maximum  chemical  absorption 
for  a given  volume. 
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a.  Shots  BB001-BB003 


b.  Shots  BB00U-BB007 


Figure  2.  Containers  used  to  hold  AN/FO 


Log  8 


9.  As  stated  in  paragraph  •*,  logs  were  used  to  represent  tree 
stems.  It  was  believed  that  logs  voul'l  simulate  stem  behavior  under 
shock  loading,  since  the  previous  study1  had  shown  that  the  failure  of 
a tree  stem  is  not  influenced  by  the  crown  or  root  mass  of  the  tree, 
i.e.,  the  portion  of  the  stem  that  is  subjected  to  deformation  is  a 
free  body  from  the  initial  interaction  of  the  stem  and  the  shock  front 
until  stem  failure.  Logs  used  for  the  test  program  were  termed  "free 
flying"  to  indicate  that  their  motion  was  in  no  way  constrained,  i.e. 
a free  body.  The  shock  front  that  accelerated  the  logs  into  free  flight 
was  assumed  to  have  spherical  symmetry;  therefore,  the  impulse  loading 
on  a free  body  (log)  was  the  same  in  any  plane  tangent  to  the  shock 
front.  All  the  logs  were  statically  placed  (prior  to  detonation)  in  a 
horizontal  position  such  that  the  longitudinal  axis  of  the  log  was 
tangent  to  the  shock  front  (assuming  that  the  shock  front  would  expand 
as  a sphere  centered  at  GZ).  The  logs  were  placed  horizontally  in  lieu 
of  other  orientations  because  they  could  be  positioned  easily  on  simple 
stands  with  only  minor  difficulties  in  maintaining  calibration  on  logs 
that  were  instrumented. 

10.  Logs  used  during  the  first  five  shots  (see  Table  2)  were 
initially  positioned  no  closer  than  8.8  m to  GZ  and  were  greater  than 
6 cm  in  diameter.  These  constraints  were  necessary  so  that  the  logs 
would  not  sustain  physical  damage  (such  that  the  log  would  not  be  re- 
usable, see  Figure  3)  and  the  instrumentation  would  not  malfunction 
because  of  abnormally  high  overpressures  and  Jetting  (see  paragraphs 
6U-66 ) . The  data  from  the  Fort  Benning  test  shot1  (3.6-metric-ton 
yield)  indicated  that  all  trees  6 cm  in  diameter  or  larger  that  were 
30  m or  more  from  GZ  did  not  fail.  According  to  the  cube-root  scaling 
laws  applicable  to  explosives,^  a 6-cm-diam  stem  would  not  fail  at 
8.8  m from  GZ  for  a 0.09-metric-ton  yield.  Therefore,  by  placing  the 
logs  at  a distance  of  8.8  m or  more  from  GZ  (an  overpressure  of  approxi- 
mately 40  psi)  and  by  using  logs  greater  than  6 cm  in  diameter,  it  was 
felt  that  no  serious  damage  to  the  logs  nor  the  instrumentation  mounted 
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on  some  of  them  would  result.  The  logs  were  moved  closer  to  GZ  for 
tests  BB006  and  BB007  to  determine  the  limit  of  reliability  of  the 
instrumentation. 


Figure  3.  Maximum  expansion  of  fireball  (test  shot 
BB002)  did  not  reach  logs  placed  8.8  m from  GZ 


Noninstrumented  logs 

11.  Two  noninstrumented  logs  were  used  in  shots  BB001  and  BB002 
to  obtain  a displacement-time  history  of  a flying  log  by  high-speed 
photography  (see  paragraph  55).  The  logs  were  placed  prior  to  detona- 
tion of  the  explosive  on  a stand  constructed  of  7-6-cm  channel  iron 
anchored  in  concrete  as  shown  in  the  extreme  right  of  Figure  3 (also  see 
Figure  ~).  The  logs  were  111.8  cm  long  and  had  20-  and  10-cm  diameters 
and  weighed  33,9 66  and  6,351  g,  respectively.  The  ends  were  painted  in 
alternating  black  and  white  quadrants. 

Instrumented  logs  and  stands 

12.  Fifteen  logs  were  instrumented;  most  were  used  for  more  than 
one  test  shot.  These  logs  were  assigned  numbers  f^r  test  purposes  as 
indicated  in  Table  2;  however,  log  number  1 in  one  test  was  not  neces- 
sarily the  same  log  number  1 in  another  test.  These  logs  were  111.8  cm 
long,  either  10  or  20  cm  in  diameter,  and  varied  in  mass  (including 
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instrumentation)  from  9,054  to  38,206  g.  Fourteen  of  the  logs  were 
rough  (natural  surface);  the  other  log  had  been  turned  down  in  a lathe 
so  that  its  surface  was  smooth  and  was  used  to  study  the  relative 
effect  of  surface  condition  on  the  drag  coefficient  of  the  log.  The  logs 
were  instrumented  on  both  ends  and  the  center  to  obtain  a displacement- 
time history  of  a flying  log  by  using  accelerometers  (paragraph  24) , 
velocity  gages  (paragraph  25),  and  specially  constructed  displacement- 
time equipment  (paragraph  26).  Because  the  displacement-time  equipment 
used  with  the  logs  is  unique,  a brief  description  of  the  several  types 
used  during  the  test  program  is  provided  in  paragraph  26. 

13.  Two  types  of  stands  were  used  during  the  test  program  to  sup- 
port the  instrumented  logs  prior  to  detonation  of  the  explosive.  The 
first  type  was  constructed  of  7.6-cm  channel  iron  anchored  in  buried 
concrete  bases  of  relatively  large  mass  (Figure  4a).  The  backfill 
around  each  base  was  compacted  to  a high  density.  This  construction 
technique  was  used  to  minimize  the  effect  of  any  vibrations  caused  by 
ground  motion  occurring  at  the  test  stands  before  the  arrival  of  the 
shock  front  and  to  prevent  any  damage  to  the  stands  by  airblast.  Exami- 
nation of  the  log  stands  and  test  data  after  the  first  test  shot  indi- 
cated that  the  rigid  construction  techniques  described  above  were  not 
necessary.  Therefore,  a less  costly  mounting  stand  was  prefabricated 
from  1.25-m  lengths  of  5.1-cm  angle  iron  (vertical  members)  and  15-cm 
sections  of  7.6-cm  channel  iron  (horizontal  members)  as  shown  in 
Figure  4b.  One  vertical  and  one  horizontal  member  were  required  to 
support  each  end.  This  stand  could  be  installed  by  one  man  whereas  the 
concrete-based  stands  required  several  men  and  a dragline  to  lift  the 
concrete  base  into  place.  All  logs  were  positioned  on  the  stands  such 
that  their  initial  trajectories  would  not  be  affected  by  electrical 
cables  or  by  the  supports  on  which  they  rested.  The  logs  were  put  in 
place  by  resting  the  ends  on  foam  pads.  Each  log  instrumented  with 
accelerometers  and/or  velocity  gages  was  staked  to  prevent  any  movement. 
The  stability  was  necessary  so  that  the  accelerometers  and  velocity 
gages  would  remain  calibrated  during  the  time  interval  between 
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a.  Stand  constructed  of  7.6-cm  angle 
iron  anchored  in  concrete 


b.  Stand  constructed  of  a 1.25-m  length  of  5.1-cm  angle  iron 
(vertical  member)  and  a 15-cm  section  of  7.6-cm  channel  iron 

(horizontal  member) 

Figure  1*.  Log  stands  used  in  test  program 
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positioning  and  detonation.  The  stakes  were  removed  Just  prior  to 
countdown  and  detonation. 

14.  An  attempt  was  made  in  shot  BB003  (paragraph  36)  to  evaluate 
the  effect  of  the  supporting  stand  on  the  test  results  of  an  instru- 
mented log;  this  was  accomplished  by  hanging  log  number  3 (see  Table  2) 
by  its  ends  with  piano  wire.  The  piano  wire  was  hung  across  two  elec- 
trodes (mounted  in  the  horizontal  beam  as  shown  in  Figure  5a)  that  were 
part  of  a circuit  connected  directly  to  the  firing  mechanism.  When  the 
charge  was  detonated,  a high  current  passed  through  the  piano  wire,  in- 
stantaneously melting  the  wire  and  releasing  the  log  into  free  fall.  The 
entire  support  array  was  constructed  of  10-  by  10-cm  timber  (Figure  5b). 

Recorders 


15.  The  field  data  (other  than  the  high-speed  photographic  rec- 
ords) from  all  the  shots  (Table  l)  were  recorded  on  magnetic  tape. 
Overpressure,  acceleration,  velocity,  displacement- time,  TOA,  a zero- 
time indication,  and  a time  reference  code,  i.e.  an  Inter-Range  Instru- 
mentation Group  (IRIG)  trace,  were  recorded  during  each  detonation. 

The  data  were  recorded  with  Model  3560  Sangamo  Portable  Magnetic  Tape 
Recorder /Reproducers  manufactured  by  Sangamo  Electric  Company.  Two 
such  recorders  were  used  for  each  detonation;  each  carried  14  data 
channels.  The  IRIG  trace  was  obtained  with  a Model  SP-100-BS  IRIG-B 
Test  Time  Code  Generator  manufactured  by  Flow  Corporation  and  fed  into 
the  recorder.  The  generator  produced  100  pps,  30  bit  binary  coded 
decimal  time  of  year/17  bit  time  of  day.  The  recorders  and  associated 
equipment  were  housed  in  a van  located  200  ra  from  GZ. 

Instrumentation 


l6.  Instrumentation  was  selected  for  each  detonation  that  would 
provide  information  on  the  blast  characteristics  and  the  response  to 
the  blast  of  free-flying  logs  (see  paragraph  9 for  definition  of  free- 
flying  logs).  The  characteristics  of  the  blast  were  determined  from 
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b.  Log  suspended  by  piano  wire  hung  across  two 
electrodes  mounted  in  overhead  beam 

Figure  5-  Apparatus  used  as  alternative  to  log  stands 
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measurements  of  incident  overpressure  versus  time,  time  of  arrival  of 
the  shock  front  at  various  distances  from  GZ,  and  from  high-speed- 
photography  records.  The  flying  log’s  response  to  the  shock  front  was 
evaluated  from  measurements  of  acceleration,  velocity,  and  displacement- 
time, and  from  high-speed-photography  records. 

Overpressure  (OP)  gages 

17.  Norwood  Models  211-3^-C  and  llA-3^  bonded  strain-gage  pres- 
sure transducers  manufactured  by  American  Standard,  Inc.,  were  used  to 
detect  the  TOA  of  the  shock  front  and  to  measure  the  OP  at  selected 
locations  (see  Figure  1 and  Table  3).  These  transducers  were  mounted 
flush  with  the  ground  (Figure  6a).  Pressure  imposed  upon  the  flush- 
mounted  diaphragm  and  strain-tube  assembly  causes  minute  compression  of 
the  strain  tube,  resulting  in  a signal  output  from  the  gage  that  is 
directly  proportional  to  the  pressure.  Also  used  were  Model  LC-35 
pressure  transducers  manufactured  by  Atlantic  Research  Corporation. 

These  gages  were  used  in  conjunction  with  a Kistler  preamplifier  to 
achieve  the  necessary  signal  conditioning  and  were  mounted  on  stands 
such  that  the  points  were  oriented  radially  toward  GZ  (Figure  6b). 

TOA  gages 

18.  The  TOA  gages,  constructed  at  WES,  were  used  in  shots  BB001, 
BB002,  and  BB005  to  determine  the  arrival  time  of  the  shock  front  at  a 
given  distance  from  GZ  (Figure  1 and  Table  h).  Each  gage  consisted  of 
a small  paper  vane  attached  to  an  arm  that  pivoted  back  against  a brit- 
tle pencil  lead  (Figure  7a).  As  the  shock  front  passed,  the  lead  was 
broken  by  the  arm  and  a voltage  change  was  recorded  on  magnetic  tape. 

The  TOA  gages  were  assembled  in  arrays  to  allow  study  of  the  symmetry 
of  the  shock  front  and  to  obtain  continuous  records  of  OP  versus 
distance  from  GZ.  Figure  7b  shows  the  array  used  in  shot  BB002  (see 
paragraph  35)  to  obtain  data  from  which  OP  could  be  calculated  from  TOA 
data  (see  paragraph  80)  and  from  which  the  shape  of  the  shock  front 
could  be  examined  (paragraph  72). 

19.  The  TOA  gages  were  mounted  such  that  the  pressure-sensitive 
element  (the  paper  vane)  was  between  GZ  and  the  aluminum  pole  on  which 
the  gage  was  mounted.  Hence,  the  pencil  lead  would  be  broken  by  the 
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b.  Atlantic  Research  Corporation 
Model  LC-35  gage  oriented  radially 
toward  GZ 

Figure  6.  OP  gages 
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a.  Single  TOA  gage 


b.  Array  of  TOA  gages  used  for  shot  BB002 
Figure  7»  TOA  gages 
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paper  vane  before  the  shock  front  set  the  aluminum  pole  into  oscillatory 
motion  (and,  as  a result,  possibly  break  the  pencil  lead).  The  TOA  gages 
were  attached  to  the  aluminum  poles  with  hose  clamps;  a Styrofoam  strip 
was  placed  between  the  hose  clamp  and  the  pole  to  mechanically  decouple 
the  gage  and  the  aluminum  pole.  This  was  necessary  because  wind  or  per- 
sonnel occasionally  set  an  aluminum  pole  into  oscillation  and  the  pencil 
lead  was  broken.  The  Styrofoam  strip  prevented  breakage  of  the  pencil 
lead  in  the  majority  of  cases  where  unwanted  oscillations  occurred. 

20.  Two  electrical  circuits  were  developed  to  accommodate  the  TOA 

gages.  One  circuit  (Figure  8a)  was  used  for  a single  gage  on  a single 
data  channel  and  the  other  (Figure  8b)  for  multiple  gages  on  a single 
channel.  The  circuit  shown  in  Figure  8a  functions  as  follows:  Prior 

to  the  breaking  of  the  pencil  lead  by  the  paper  vane,  the  amplifier 
draws  no  current  from  the  2U-v  source.  When  the  pencil  lead  breaks, 
the  input  voltage  to  the  amplifier  rises  <o  5.1  v (the  operating  volt- 
age of  the  Zener  diode)  and  the  amplifier  draws  current.  The  increase 
in  current  produces  a discontinuity  on  a magnetic  tape  data  record. 

This  discontinuity  can  be  visually  examined  on  an  oscillogram  reproduced 
fiom  the  magnetic  tape.  The  arrival  time  can  be  estimated  from  the 
discontinuity.  In  addition  to  its  capacity  as  a voltage  regulator,  the 
Zener  diode  serves  as  a clipper  for  transients  induced  in  the  cable  by 
the  explosion.  This  circuit  was  used  when  several  TOA  values  were 
required  at  the  same  distance  from  GZ,  i.e.,  separate  channels  were 
required  because  the  discontinuities  occurred  simultaneously. 

21.  The  circuit  shown  in  Figure  8b  functions  as  follows:  Prior 

to  the  breaking  of  the  pencil  lead  closest  to  GZ,  the  input  load  imped- 
ance on  the  amplifier  is  the  equivalent  resistance  of  the  number  of 
10— k resistors  in  parallel  (one  10-k  resistor  for  each  TOA  gage).  The 
input  voltage  is  zero.  When  the  shock  front  breaks  the  first  pencil 
lead,  the  Zener  diode  becomes  active  and  supplies  5.1  v to  the  10-k 
resistor  in  series  with  the  amplifier  input,  and  the  amplifier  draws 
current,  which  in  turn  produces  a discontinuity  on  a magnetic  tape  data 
record  that  corresponds  to  the  arrival  time  of  the  shock  front.  When 
the  second  pencil  lead  breaks,  the  amplifier  draws  twice  as  much  current. 


21 


PENC'L 

LEAD 


TO  AMPLIFIER 


a.  Circuit  for  one  TOA  gage  on  single  data  channel 
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b.  Circuit  for  multiple  TOA  gages  on  single  data  channel.  Only- 
three  TOA  gages  are  shown,  but  the  circuit  can  be  expanded  to 
accommodate  more  gages  (see  text) 


Figure  8.  Circuits  for  TOA  gages 


and  another  discontinuity  occurs,  i.e.,  the  currents  from  the  individ- 
ual TOA  circuits  are  additive.  The  increase  of  current  due  to  the 
breaking  of  the  pencil  leads  continues  until  the  shock  front  reaches  the 
end  of  the  array  of  TOA  gages.  The  number  of  gages  on  a circuit  is  lim- 
ited only  by  the  readability  of  the  record,  i.e.,  as  the  number  of 
gages  on  the  line  increases,  the  step  height  of  the  discontinuities 
decreases  for  the  same  maximum  amplifier  input  current.  The  maximum 
number  of  gages  successfully  used  on  a single  data  channel  during  the 
test  program  was  15.  This  type  of  circuit  was  used  when  TOA  data  were 
required  on  the  same  azimuth  as,  but  at  different  distances  from,  GZ. 

Instrumentation  for 
determination  of  zero  time 

22.  The  estimation  of  shock  front  arrival  time  from  oscillogram 
trace  discontinuities  requires  an  accurate  reference  value  (zero  time). 
Zero  time  was  generated  by  two  techniques.  The  primary  zero  time  indi- 
cation was  obtained  from  a break-wire  circuit.  The  wire  was  wrapped 
around  the  booster  in  the  center  of  the  explosive  sphere  of  AN/FO  (para- 
graph 8).  When  the  explosive  was  detonated,  the  wire  broke  and  produced 
a change  of  voltage  in  the  circuit.  The  change  of  voltage  appeared  as 

a sharp  discontinuity  on  the  resulting  oscillogram  trace.  The  secondary 
zero  time  was  obtained  from  a Questar  telescope  (see  Table  l)  with  a 
photodiode  attached  to  the  eyepiece  (Figure  9).  The  illumination  from 
the  explosion  generated  sufficient  output  from  the  photodiode  to  place 
a discontinuity  on  the  oscillogram  trace. 

Earth-movement  canister 

23.  Horizontal  and  vertical  soil  displacements  were  measured  by 
a Model  35^137-0102  velocity  gage  manufactured  by  Bell  and  Howell.  The 
gage  was  used  only  for  test  shots  BB001  and  BB002  (see  Table  l).  Data 
from  these  tests  indicated  that  the  air  shock  was  arriving  before  the 
ground  shock  and,  therefore,  ground  shock  could  not  affect  the  instru- 
mentation (see  paragraph  6l);  hence,  use  of  these  gages  was  discon- 
tinued after  test  shot  BB002. 

Accelerometers 

2k.  Models  226l  (hOOO-g  max)  and  2265M1  (2500-g  max) 
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Figure  9-  Questar  telescope  with  photodiode  attached 
to  eyepiece  used  to  obtain  secondary  zero  time 


accelerometers  manufactured  by  Endevco  Corporation  were  used  (see 
Table  2)  to  measure  the  horizontal  and  vertical  components  of  the  accel- 
eration experienced  by  the  flying  logs.  The  choice  of  accelerometer 
range  is  usually  made  on  the  basis  of  the  maximum  acceleration  expected. 
However,  the  range  choice  also  dictates  the  resonant  frequency  of  the 
accelerometer  (a  high-range  accelerometer  has  higher  resonant  frequency 
than  a low-range  accelerometer  of  the  same  type).  The  Model  226l  accel- 
erometer had  a resonant  frequency  of  30  kHz  and  the  Model  226UM1,  20  kHz. 
To  avoid  getting  erroneous  acceleration  data,  the  gages  must  be  operated 
well  below  their  resonant  frequencies,  as  was  the  case  for  this  test 
program.  The  accelerometers  were  mounted  within  hollowed-out  portions 
of  the  flying  logs,  and  expanded  foam  was  used  to  secure  each  gage 
in  place  (Figure  10).  The  accelerometers  produced  signals  with  ampli- 
tudes proportional  to  acceleration  along  the  designated  axis  (horizon- 
tal or  vertical)  of  the  log.  They  were  used  in  conjunction  with  a 
WES  !fyPe  100-0P  amplifier  to  achieve  the  necessary  signed  conditioning 
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Figure  10.  Acceleration  and  velocity  gages  mounted 
in  expanded  foam  in  holloved-out  portion  of  log. 

Cables  associated  with  the  gages  are  shown  coming 

from  the  log 

and  were  calibrated  prior  to  each  test  shot. 

Velocity  gages 

25*  Two  types  of  velocity  gages  were  used:  a Model  7LV6-N  veloc- 

ity gage  manufactured  by  Sanborn  Company,  and  a Model  DX5  velocity 
gage  manufactured  by  Sparton  Southwest,  Inc.  The  Sanborn  gage  (a  linear 
inductive  type)  produced  a signal  with  an  amplitude  proportional  to 
linear  velocity  and  had  an  effective  measuring  length  (displacement)  of 
15  cm.  This  gage  was  used  on  log  numbers  2 and  3 in  test  shots  BB003 
and  BB004  (see  Table  2).  The  Sparton  gage  (variable  reluctance  type) 
produced  a signal  with  an  amplitude  proportional  to  horizontal  or  verti- 
cal velocity.  This  gage  was  used  on  log  numbers  1 and  3 in  test  shots 
BB003  and  BBOOU.  Each  gage  was  precalibrated  (or  checked  for  calibra- 
tion) prior  to  each  test  shot.  The  necessary  temperature  correction 
factors  (paragraph  29b(6))  were  incorporated  in  the  calibration  proce- 
dure. The  velocity  gages  were  mounted  in  the  logs  in  the  same  manner 
as  the  acceleration  gages  (Figure  10). 

Displacement-time  gages 

26.  Four  types  of  gages  were  constructed  and  tested  at  WES  to 
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examine  the  displacement-time  characteristics  of  flying  logs.  Because 
these  gages  are  unique,  their  development  is  discussed  below.  The  ini- 
tial configuration  consisted  of  a rigid  steel  wire  connected  to  the  end 
of  the  log.  The  other  end  of  the  wire  was  hook-shaped;  perpendicular 
to  the  hook  were  fine  wires  spaced  a known  distance  apart  (Figure  11a). 
When  the  log  was  acted  upon  by  the  shock  front,  it  was  displaced.  As 
the  log  was  displaced,  the  hook  broke  the  fine  wires  in  succession. 

Figure  lib  is  a schematic  of  the  circuit.  Each  broken  wire  increased 
the  input  resistance  (and  voltage)  of  the  data  channel  amplifier.  This 
produced  a step  on  the  resulting  oscillogram  trace  corresponding  to  a 
known  distance.  Hence,  in  conjunction  with  the  IRIG  trace  (see  para- 
graph 15),  a displacement-time  record  was  obtained.  Unfortunately,  the 
apparatus  generated  data  of  poor  quality.  Even  though  the  fine-wire  con- 
figuration was  mounted  under  protective  cover  in  a rigidly  mounted,  stand 
(see  paragraph  13),  some  of  the  fine  wires  were  broken  by  the  initial 
impact  of  the  shock  front  before  the  steel  hook  came  into  contact  with 
them.  Figure  Ua  shows  the  log  in  position  on  the  stand;  the  fine-wire 
configuration  is  inside  the  hollow  upper  horizontal  member  of  the  stand. 
This  device  was  used  for  test  shots  BB001  and  BB002  (see  Table  2).  The 
fine-wire  concept  was  rejected  in  favor  of  using  photoelectric  cells 
in  place  of  the  fine  wires.  The  circuitry  required  for  the  photo- 
electric scheme  became  extremely  complicated  and  was  subsequently  re- 
jected. The  third  configuration  consisted  of  a wire  magnetized  at  known 
increments,  attached  to  the  end  of  the  log.  As  the  log  moved  away  from 
the  stand,  the  wire  moved  past  a recording  head  mounted  in  the  stand. 
Thus,  a displacement-time  record  was  produced  for  the  flying  log.  This 
configuration  was  successful,  as  discussed  in  the  following  paragraph. 

27.  The  magnetic  recording  head  and  magnetized  wire  were  mounted 
in  the  upper  horizontal  member  of  the  stand  (Figure  l*a)  in  much  the 
same  manner  as  the  wire-breaker  and  photoelectric  configurations.  Fig- 
ure 12  is  a defiled  diagram  of  the  magnetic  head  and  magnetized  wire. 

One  end  of  the  magnetized  wire  is  attached  to  the  object  to  be  dis- 
placed and  the  other  end  is  threaded  through  the  wire  guide  on  the  mag- 
netic head  (the  head  remains  stationary).  When  the  wire  is  moved 
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a.  Schematic  of  apparatus 
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CORRESPONDING  / ▼ 

DISPLACEMENT 

b.  Schematic  of  circuit 


Figure  11.  Wire-breaker  apparatus  used  to  obtain  displacement- 
time record  of  flying  log 

Xl 


OSCILLOGRAM  OF 


RECORDING  HEAD 

Figure  12.  Magnetic  apparatus  for  obtaining  displacement 
measurements  of  log  motion 


(i.e.,  log  displacement),  the  magnetized  spots  travel  across  the  air  gap 
on  the  magnetic  head  causing  a varying  flux,  $ , through  the  core.  An 
electric  potential  is  produced  in  the  coil  of  wire  proportional  to  the 
rate  of  change  in  flux,  i.e.,  V = d<J>/dt  , and  this  voltage  is  fed  into 
an  amplifier  to  raise  the  signal  to  recording  level.  A displacement- 
time trace  can  then  be  reproduced  on  an  oscillogram  for  data  analysis. 
This  device  was  used  for  test  shots  BB003  and  BBOOI4  (see  Table  2). 

28.  The  fourth  displacement-time  apparatus  is  shown  in  Figure  13. 
The  log  was  positioned  on  the  test  stand  constructed  from  1.25-m  lengths 
of  angle  iron  (vertical  members)  and  15-cm  sections  of  channel  iron 
(horizontal,  members)  (see  paragraph  13  and  Figure  Vb).  The  magnetized- 
wire  displacement- time  apparatus  (paragraph  27)  was  mounted  inside  a 
6.3-cm-diam  galvanized  pipe  structure  immediately  behind  the  test  stand. 
A protective  cap  was  threaded  over  the  open  end  of  the  pipe  to  protect 
the  magnetic  head  from  the  airblast.  The  magnetized  wire  was  threaded 
through  a Teflon  guide  in  a hole  drilled  in  the  protective  cap  (Fig- 
ure lU).  This  configuration  produced  good  quality  displacement-time 
data  and  was  easier  to  install  than  any  of  the  other  methods  previously 
tried,  and  was  used  for  test  shots  BB003  through  BB007. 

Error  analysis 

29.  Errors  are  inherent  in  any  measurement  system.  The  errors 
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Figure  13.  Modified  stand  used  with  displacement-time 
apparatus  in  test  shots  BB003-BB007.  Vertical  boards 
in  background  are  for  studying  the  shape  of  the  shock 
front  using  high-speed  photography 


introduced  by  the  gages  used  in  the  test  series,  together  with  esti- 
mates of  error,  are  as  follows: 

a_.  Acceleration  gages 

( 1 ) Calibration  error  (steady-state  response).  The 
acceleration  gages  were  calibrated  on  a rotating 
table  prior  to  their  use  in  the  field.  Each  had 
its  own  calibration  factor;  unless  the  gage  is 
operated  beyond  its  maximum  range,  this  factor  re- 
mains approximately  the  same  throughout  its  life. 
Anticipated  error  is  +3  percent. 

(2)  Frequency  response  error.  The  acceleration  gages 
have  a damping  factor  of  0.03  (undamped)  and 
therefore  have  a high  peak  in  the  response  curve  at 
the  resonant  frequency.  To  avoid  erroneous  data, 
the  gage  must  be  operated  well  below  the  resonant 
frequency.  Large  errors  can  also  occur  if  the 
coupling  device  or  medium  produces  secondary  peaks 
in  the  response  curve.  Anticipated  error  is  jOLO 
percent. 

(3)  Stability  error.  An  acceleration  gage  may  show  a 
change  in  the  calibration  factor  from  day  to  day. 
This  is  a long-term  drift  factor  that  depends  pri- 
marily on  temperature.  Anticipated  error  is 

+1  percent. 
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(4)  Field  reference  error.  The  full-scale  range  of  an 
acceleration  gage  is  obtained  from  amplifier  cali- 
bration by  a resistor  substitution  technique.  This 
resistance  can  be  matched  only  approximately.  Antic- 
ipated error  is  +1  percent. 

(5)  Linearity  error.  The  calibration  curves  are  assumed 
to  be  linear,  though  in  practice  they  have  some 
second-order  variations.  Using  a high-range  accel- 
eration gage  for  a lov-range  recording  exaggerates 
this  error.  Anticipated  error  is  +3  percent. 

(6)  Mounting  error.  The  acceleration  gages  were  "soft" 
mounted  in  the  logs,  as  opposed  to  "hard"  mounting 
in  which  the  gages  are  attached  directly  to  the  log. 
The  soft  mount  requires  that  the  transducer  be  em- 
bedded in  a lossy  material  (foam  plastic  or  poly- 
styrene) in  a hole  in  the  log.  This  allows  only  the 
low-frequency  components  of  the  log's  acceleration 
to  couple  to  the  acceleration  gage.  In  effect,  the 
soft  mount  is  a low-pass  filter.  Thus,  the  acceler- 
ation curve  or  the  integrals  of  the  curve  do  not 
display  the  log  motion,  but  represent  the  motion  of 
the  accelerometer.  Errors  in  the  acceleration  mag- 
nitude may  be  due  also  to  the  gage  not  being  exactly 
in  the  vertical  or  horizontal  plane.  Amplitude  re- 
duction is  proportional  to  the  cosine  of  the  angle 
(5-  to  10-deg  max).  Anticipated  error  is  +1  percent. 

(7)  Recording  and  digitization  error.  Minimal  error  is 
introduced  by  the  analog  tape  recorder  and  the  ana- 
log to  digital  converter  (see  paragraph  50).  Antic- 
ipated total  error  is  +0.25  percent. 

(8)  Integration  error.  If  errors  in  the  measuring 
equipment  cause  the  vecordings  of  acceleration  to 
show  an  offset  after  the  main  or  peak  acceleration 
has  been  passed,  the  integral  of  these  data  (especi- 
ally for  long  periods  after  the  peak  acceleration) 
will  contain  errors  in  proportion  to  the  product  of 
the  offset  and  the  integration  time.  The  offset 
voltage  can  be  due  to  hysteresis  in  the  amplifier 

or  in  the  acceleration  gage.  The  anticipated  error 
due  to  hysteresis  in  the  amplifier  is  j_0.1  percent, 
and  in  the  acceleration  gage  is  +1  percent. 

b.  Velocity  gages.  Comments  on  the  acceleration  gages  in 
subparagraphs  (1 ) , ( 4 ) , ( 5 ) , and  ( 7 ) above  apply  to  the 
velocity  gages  for  error  comments  (l),  (2),  (3),  and  (4), 
respectively,  below. 

(1)  Calibration  error.  Anticipated  error  is  +2  percent. 

(2)  Field  reference  error.  Anticipated  error  is  +1 
percent. 
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(3)  Linearity  error.  Anticipated  error  is  +3  percent. 

(!+)  Recording  error.  Anticipated  error  is  +0.25  percent. 

(5)  Frequency  response  error.  The  velocity  gages  are 
highly  damped  (mechanically  damped  with  a viscous 
fluid)  and  have  a frequency  response  of  from  3 to 
300  Hz.  Since  they  are  so  highly  damped,  they  do 
not  normally  have  a peak  in  their  response  curve 
unless  it  is  produced  in  the  external  mounting. 
Anticipated  error  is  +15  percent. 

(6)  Stability  error.  Velocity  gages  are  sensitive  to 
temperature  and  must  be  corrected  for  changes  from 
the  calibration  temperature  (l  percent/°F).  With 
temperature  correction,  the  maximum  anticipated 
error  is  +2  percent. 

(7)  Mounting  error.  Velocity  gages  are  very  sensitive 
to  tilt  (a  6-deg  change  is  a full-scale  effect). 

These  gages  were  zeroed  prior  to  the  detonation  by 
carefully  changing  the  log  position  for  zero 
response.  If  any  log  rotation  takes  place  (changes 
of  1 deg  or  more),  the  effect,  on  the  output  is  sig- 
nificant. Rotation  is  evident  on  the  recording  as 
a steadily  increasing  or  decreasing  amplitude  after 
the  initial  acceleration.  Anticipated  error  is 

+1  percent. 

£.  Displacement-time  gages.  Comments  on  the  acceleration 

gages  in  subparagraphs  a(l)  and  a(7)  above  apply  to  the 

displacement-time  gages  for  error  comments  (l)  and  (2), 

respectively,  below. 

(1)  Calibration  error.  Anticipated  error  is  +2  percent. 

(2)  Recording  error.  Anticipated  error  is  +0.25  percent. 

(3)  Mounting  error.  The  displacement-time  gages  were 
mounted  in  a fixed  stand  with  one  end  of  the  record- 
ing wire  fixed  to  the  log.  A pull  in  My  direction 
(i.e.  combined  vertical  and  horizontal  displacement) 
activates  the  gage.  Also,  an  air  pressure  differen- 
tial between  the  inside  of  mount  and  the  surrounding 
area  may  force  the  recording  wire  past  the  recording 
head  without  the  corresponding  log  movement.  Antic- 
ipated error  is  +1  percent. 

Meteorological  instrumentation 

30.  Temperature  and  barometric  pressure  were  measured  prior  to 
each  test  shot.  These  data  are  summarized  in  Table  5.  The  temperature 
was  measured  with  a Friez  hygrothermograph,  Model  59**;  barometric  pres- 
sure measurements  were  obtained  from  the  weather  station  at  WES.  Wind 
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velocity  was  zero  for  each  test.  Shots  were  conducted  on  days  when 
lighting  was  optimum  for  high-speed  photography.  This  requirement  pre- 
cluded problems  associated  with  detonation  of  explosives  under  low  cloud 
ceilings. 


High-Speed  Camera 

31.  Five  high-speed  l6-mm  movie  cameras  were  used  during  the 
test  series  (see  Table  1):  a Milliken  Model  DBM  SC  manufactured  by  the 

D.  B.  Milliken  Company  (framing  rate:  1+00  frame s/ sec ond ) ; a Reflex 

manufactured  by  Kodak  (6U  f/s);  a Nova  Model  16-3  manufactured  by  Photo- 
Kinetics,  Inc.  (framing  rate  variable  to  1+0,000  f/s);  a Model  HSl+01  man- 
ufactured by  Fairchild  Camera  and  Instrument  Company  (2000  f/s);  and  a 
Fastex  Model  WFl*T  manufactured  by  Red  Lake  Laboratories,  Inc.  (5000  f/s). 
These  cameras  were  used  to  obtain  the  displacement-time  history  of  fly- 
ing logs  and  to  determine  the  shape  of  the  shock  front.  Each  camera 
was  equipped  with  dual  timing  lenses  and  gave  full-frame  rather  than 
split-frame  type  of  coverage.  The  exposure  rates  used  over  the  series 
of  shots  ranged  from  6h  to  8000  f/s.  All  cameras  were  mounted  on  tri- 
pods and  were  protected  from  flying  debris  with  sandbags.  Each  tripod 
was  rigidly  fixed  by  a vertical  tension  cable  between  the  tripod  head 
and  a ground  anchor.  Cameras  were  located  approximately  50  m from  GZ. 
Electrical  timing  circuits  automatically  activated  the  cameras  prior  to 
detonation  allowing  sufficient  time  for  the  film  to  accelerate  to  the 
required  framing  rate. 


Test  Program 


Test  shots 

32.  The  objectives  and  site  layout  are  described  for  each  test 
shot  in  the  following  paragraphs. 

33.  Test  Shot  BB001.  The  objectives  of  shot  BB001  were  to: 

a.  Obtain  the  displacement-time  history  for  flying  logs. 

b.  Examine  the  geometrical  properties  of  the  shock  front. 
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£.  Make  OP  measurements  to  verify  that  the  desired  yield  was 
obtained. 

d.  Measure  soil  displacement  induced  by  the  shock  front. 

31*.  The  displacement-time  histories  of  the  flying  logs  were  ob- 
tained with  the  wire-breaker  apparatus  described  in  paragraph  26  and 
from  vigh-speed  photography  (Table  l).  Two  logs  initially  positioned 
at  8.8  m from  GZ  were  used  for  this  test  (Table  4).  The  geometry  of 
the  shock  front  was  examined  in  the  horizontal  plane  by  placing  four 
TGA  gages  at  8.8  m from  GZ,  but  spaced  90  deg  apart.  Equal  arrival 
times  at  the  same  distances  from  GZ  would  indicate  that  the  shock  front 
was  propagating  away  from  GZ  at  the  same  velocity  in  four  directions. 

The  geometry  of  the  shock  front  was  examined  in  the  vertical  plane  using 
high-speed  photography  and  also  by  placing  three  TOA  gages  at  10.8,  12.3, 
and  13.8  m from  GZ  at  the  same  elevation  as  that  of  the  center  of  the 
explosion;  three  other  gages  were  placed  2.37  m above  these  three.  The 
inclination  of  the  shock  front  with  the  horizontal  could  be  examined 
from  the  activation  sequence  of  the  TOA  gages  (see  paragraph  71).  Four 
OP  gages  were  placed  at  8.8  m from  GZ  at  azimuths*  of  23,  90,  l8o,  and 
270  deg  to  verify  that  the  desired  yield  had  been  obtained  (Table  3).  An 
earth-movement  canister  was  placed  at  8.8  m from  GZ  to  measure  horiz- 
ontal and  vertical  soil  displacement.  Zero  time  was  obtained  from  a 
break-wire  device  and  a Questar  telescope  as  described  in  paragraph  22. 
These  two  devices  were  used  during  the  entire  test  series.  Meteorolog- 
ical data  were  taken  as  described  in  paragraph  30.  These  data  were 
taken  for  each  test  shot.  The  explosive  was  placed  in  a cubical  Styro- 
foam container  at  GZ.  A wide  variety  of  cameras  was  used  during  the 
test  series  (see  Table  l),  and  the  selection  of  a particular  high-speed 
camera  for  a particular  shot  depended  more  on  its  availability  than  on 
its  capabilities.  Some  deterioration  in  the  quality  of  the  high-speed 
photography  due  to  smoke  from  the  explosion  and  ground  dust  was  noted 
in  this  shot  (see  Figure  15). 

* Azimuth  was  relative  to  point  in  center  of  large  concrete  block  that 
was  base  of  wire-breaker  displacement-time  stand  used  in  test  shots 
BB001  and  BB002.  The  origin  of  the  coordinate  system  was  at  GZ. 


34 


Figure  15.  Ground  dust  beginning  to  obscure 
motion  of  log  (shot  BBOOl) 

35.  Test  shot  BB002.  The  objectives  of  shot  BB002  were  the  same 
as  those  of  BBOOl.  The  wire-breaker  apparatus  was  again  used  to  deter- 
mine the  displacement-time  history  of  two  flying  logs  initially  posi- 
tioned at  8.8  m from  GZ.  In  addition  to  the  four  TOA  gages  positioned 
at  8.8  m from  GZ  (spaced  90  deg  apart)  that,  were  used  to  study  shock- 
front  geometry  in  the  horizontal  plane,  four  additional  TOA  gages  were 
placed  at  22.1*8  m from  GZ  and  spaced  90  deg  apart.  Two  linear  strings 
of  eighteen  TOA  gages  were  placed  radially  from  GZ  starting  at  1.65  m 
and  extending  out  to  21. 64  m.  The  strings  were  positioned  90  deg  apart 
(see  Table  1*).  The  arrival  times  measured  at  these  gages  could  be  used 
to  compute  OP  as  a function  of  distance  from  GZ  (see  paragraph  80). 

The  arrival  times  could  also  be  used  with  the  four  gages  at  8.8  m and 
at  22.1*8  m to  study  shock  front  geometry  in  the  horizontal  plane  be- 
cause each  string  had  a gage  at  the  same  distance  from  GZ.  An  addi- 
tional string  of  TOA  gages  was  placed  1 ra  above  one  of  the  strings 
described  above  beginning  at  8.0  m from  GZ  and  extending  out  to  16.1*  m. 
This  configuration  was  used  to  examine  the  inclination  of  the  shock 
front  with  the  horizontal.  Two  OP  gages  were  placed  at  5.1*9  m from  GZ 
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(azimuths  of  93  and  183  deg,  respectively)  and  two  were  placed  at  8.8  m 
(azimuths  of  23  and  273  deg,  respectively)  (Table  3).  An  attempt  was 
made  to  reduce  ground  dust  by  placing  a membrane  over  the  ground  at  dis- 
tances close  to  GZ  that  were  in  the  field  of  view  of  the  high-speed 
cameras  (Figure  16).  The  membrane  proved  ineffective  and  was  not  used 
in  any  of  the  following  tests. 


Figure  16.  Membrane  used  for  shot  BB002  to  reduce  ground  dust 


36.  Test  shot  BB003.  The  displacement-time  instrumentation  used 
for  shots  BB001  and  BB002  proved  to  be  unreliable;  therefore,  the  major 
objective  of  shot  BB003  was  to  obtain  reliable  displacement-time  data 
for  the  flying  logs.  Three  logs  were  used  in  shot  BB003.  These  logs 
were  positioned  at  8.8  m from  GZ.  Accelerometers  and  velocity  gages 
were  mounted  on  each  log.  The  displacement-time  instrumentation  used 
in  the  first  two  shots  was  replaced  by  the  magnet i zed-wire  devices  de- 
scribed in  paragraphs  26  and  27.  One  of  the  magnet i zed-wire  devices 
was  mounted  in  the  stand  used  for  the  wire  breaker;  the  other  was 
mounted  in  a less  rigid  stand  constructed  of  6.3-cm-diam  galvanized 
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pipe  (see  paragraph  28).  No  TOA  gages  were  used  for  this  test  shot. 
Three  OP  gages  were  placed  at  5>^9  m from  OZ  at  azimuths  of  93,  183, 
and  273  deg  (Table  3).  The  technique  of  hanging  a log  by  piano  wire 
(log  number  3;  see  paragraph  14)  resulted  in  no  detectable  improvement 
in  the  quality  of  the  data  and  so  was  deleted  from  the  program.  The 
cost  and  effort  involved  in  the  installation  could  not  be  Justified. 

37.  Test  shot  BBOO^.  The  objective  of  this  shot  was  to  collect 
additional  displacement-time  data  for  the  flying  logs.  Four  logs  were 
positioned  at  8.8  m for  this  shot.  One  log  was  attached  to  the 
magneti  zed-wire  device  mounted  in  the  stand  used  for  the  wire  breaker 
in  shots  BB001  and  BB002.  The  other  three  were  attached  to  the 
magneti  zed-wire  apparatus  mounted  in  the  6.3-cm-diam  galvanized  pipe 
stand.  Accelerometers  and  velocity  gages  were  mounted  on  three  logs. 

The  fourth  log  had  only  an  accelerometer.  This  log  had  a smooth  sur- 
face and  was  used  to  examine  the  effect  of  the  log's  surface  on  the 
displacement-time  history.  All  other  logs  used  during  the  test  series 
had  their  natural  surface  (see  paragraph  12).  No  TOA  gages  were  used 
for  this  test.  Three  OP  gages  were  placed  at  5.^9  m from  GZ  at  azimuths 
of  138,  183,  and  273  deg  (Table  3).  Reasonable  success  had  been  experi- 
enced after  shots  BB002  and  BB003  in  locating  the  position  of  the  shock 
front  at  various  times  after  detonation  by  viewing  adjacent  frames  of 
the  high-speed  photography  (see  paragraph  56).  To  improve  the  contrast 
of  the  background,  three  plywood  backboards  (1.2  by  2.U  m)  were  erected 
in  the  field  of  view  of  the  high-speed  cameras  (see  Figure  13).  These 
boards  were  positioned  perpendicularly  to  the  camera  and  were  painted 
with  black  and  white  diagonal  stripes;  they  were  used  for  the  remainder 
of  the  test  series.  The  cubical  Styrofoam  container  (Figure  2a)  used 

to  hold  the  AN/FO  prior  to  detonation  for  test  shots  BB001  through  BB003 
was  replaced  by  a spherical  container  (Figure  2b  and  Table  l)  to  pre- 
clude the  possibility  that  the  additional  volume  of  the  cubical  con- 
tainer might  contribute  to  the  heat  of  the  explosion.  The  spherical 
container  was  used  for  the  remainder  of  the  test  series. 

38.  Test  shot  BB005.  The  magnetic-wire  displacement-time  appara- 
tus mounted  in  the  6.3-qm-diam  galvanized  pipe  stand  generated  excellent 
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data  in  shot  BB004,  and  was  used  for  the  remainder  of  the  test  shots. 

Now  that  reliable  instrumentation  had  been  assembled,  a test  shot  (BB005) 
was  conducted  to  hopefully  meet  all  the  design  objectives  of  test  shot 
BB001.  Two  logs  were  initially  positioned  at  8.8  and  12.2  m,  respec- 
tively, from  GZ  (Table  2).  All  logs  in  prior  tests  had  been  positioned 
at  8.8  m;  the  log  placed  at  12.2  m was  used  to  verify  that  the  magnetic- 
wire  displacement-time  device  would  function  correctly  when  the  incident 
OP  on  the  log  was  less  than  40  psi  (this  was  the  predicted  OP  at  8.8  m 
for  the  104.5-kg  AN/FO  charge).  The  predicted  OP  at  12.2  m was  20  psi. 
Five  TOA  gages  were  placed  at  8.8  m from  GZ  at  azimuths  of  68,  203,  247, 
292,  and  3^2  deg,  respectively  (Table  4).  A linear  string  of  19  gages 
was  oriented  at  150  deg  beginning  at  4.46  m from  GZ  and  extending  out  to 
22.48  m.  An  experiment  was  conducted  during  this  test  to  assess  the 
effects  of  ground  motion  on  the  TOA  gages.  Three  TOA  gages  were  placed 
in  cylindrical  cavities  at  6.51,  9-22,  and  13.65  m from  GZ.  The  covers 
of  these  cavities  were  sealed  flush  with  the  ground.  Six  OP  gages  were 
used  for  this  shot  (Table  3),  placed  at  distances  of  5.49,  5-49,  8.8, 

8.8,  8.8,  and  12.2  m from  GZ  at  azimuths  of  108,  293,  207,  344,  331, 
and  303  deg,  respectively. 

39-  Test  shot  BBQ06.  One  log  was  used  for  this  test.  It  was  po- 
sitioned at  4.4  m from  GZ  to  verify  that  the  magnetic-wire  displacement 
device  would  function  correctly  at  an  OP  of  200  psi,  and  that  the  gal- 
vanized pipe  stand  would  not  fail.  No  TOA  gages  were  used  for  this 
test.  OP  gages  were  placed  at  distances  of  5.49,  5-49,  and  8.8  m from 
GZ  at  azimuths  of  108,  207,  and  292  deg,  respectively  (Table  3). 

40.  Test  shot  BB007.  This  test  was  identical  with  shot  BB006  ex- 
cept that  the  log  was  positioned  at  2.2  m from  GZ  (a  predicted  OP  of 
700  psi).  The  galvanized-pipe  support  apparatus  after  detonation  of 
the  charge  is  shown  in  Figure  17* 

General  procedures 

41.  Prior  to  each  shot,  all  items  that  might  have  affected  the 
trajectory  of  the  flying  logs  or  might  have  caused  malfunction  of  the 
OP  or  TOA  gages  were  removed  from  the  test  area.  This  included  burying 
all  cables,  raking  up  all  loose  grass,  paper,  etc.,  and  removing  all 
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Figure  17.  Gal vani zed-pipe  support  apparatus 


stakes  that  held  the  instrumented  logs  in  place. 

42.  Immediately  after  the  detonation,  the  entire  test  area  was 
surveyed.  This  included  surveying  the  final  positions  of  the  flying 
logs  and  all  fixed  instrumentation  apparatus.  Survey  of  the  apparatus 
was  necessary  because  much  of  the  data  required  that  the  support  appara- 
tus and  gages  be  fixed  with  respect  to  a common  frame  of  reference; 
otherwise,  the  data  could  not  be  related. 

43.  Although  not  directly  related  to  the  objectives  of  the  test 
program,  the  following  information  may  be  of  interest.  Test  shot  BB001 
produced  a crater  60  cm  deep  and  285  cm  wide  as  shown  in  Figure  l8a. 

The  dimensions  of  the  crater  indicated  that  some  of  the  explosive  yield 
was  diverted  from  the  shock  front  to  cratering.  This  problem  was  elimi- 
nated by  filling  the  crater  after  each  subsequent  explosion  with  a mate- 
rial that  was  denser  and  more  cohesive  than  the  original  soil.  This 
material  was  a soil-cement  consisting  of  10  percent  cement  by  weight. 

The  crater  was  backfilled  by  compacting  the  soil-cement  in  successive 
lifts  of  15  to  20  cm.  Use  of  this  technique  reduced  the  crater  to  9 cm 
deep  and  200  cm  wide  for  shot  BB002  (Figure  l8b). 


39 


b.  Test  shot  BB002 

Figure  18.  Craters  produced  by  test  shots  BB001  and  BB002 
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PART  III:  PROCESSING  OF  EXPERIMENTAL  DATA 


L4.  All  experimental  data  (except  high-speed  photography)  were 
recorded  on  magnetic  tape  and  reproduced  on  oscillograms  for  analysis. 
These  data  included  shock-front  TOA  and  peak  OP,  and  the  dynamic  prop- 
erties of  the  flying  logs.  Prior  to  analysis,  all  data  vere  processed 
using  the  methods  described  in  the  following  paragraphs.  A block 
diagram  shoving  the  data  processing  procedure  is  shown  in  Figure  19. 

Di scontinuouB  Data  Records 


Time  of  arrival 

h5.  TOA  data  vere  obtained  directly  from  the  oscillograms  pro- 
duced from  magnetic  tape.  Arrival  time  of  the  shock  front  at  a partic- 
ular TOA  gage  produced  a step  discontinuity  on  the  oscillogram  trace. 
This  discontinuity  resulted  from  a change  in  the  circuit  load  resistance 
caused  by  the  breaking  of  the  pencil  lead  (see  paragraphs  18-21). 

L6.  A value  for  the  arrival  time  was  obtained  by  projecting  a 
line  from  the  discontinuity  to  the  IRIG  trace  (this  line  must  be  perpen- 
dicular to  the  baseline  of  both  traces).  The  arrival  time  was  then 
estimated  by  counting  the  number  of  cycles  on  the  IRIG  trace  between 
zero  time  and  the  intersection  of  the  projected  line  from  the  disconti- 
nuity (Figure  20). 

Log  displacement-time  history 

U7.  The  flying  log  displacement-time  history  was  obtained  by 
using  the  magnetized  steel  wire  and  recording  head  apparatus  described 
in  paragraphs  26-28.  The  magnetized  spots  on  the  steel  wire  vere  2 cm 
apart;  as  a spot  passed  the  recording  head,  a voltage  spike  was  produced. 
The  resulting  oscillogram  trace  is  a series  of  spikes  (Figure  21).  To 
evaluate  the  trace,  a line  drawn  perpendicular  to  the  trace  was  pro- 
jected from  the  spike  to  the  baseline  constructed  from  the  IRIG  trace, 
thereby  relating  displacement  to  time.  By  following  this  procedure 
for  each  spike,  a time-dependent  history  of  log  motion  was  produced. 
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Figure  19.  Data  processing  procedure 
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Figure  21.  Method  used  to  relate  displacement  and 
time  using  magnet i zed-wire  apparatus  (from  test 
shot  BB006).  Time  shown  on  baseline  is  the  time 
after  detonation 

Continuous  Data  Records 

U8.  Peak  OP  and  the  dynamic  parameters  associated  with  the  motion 
of  the  instrumented  logs  were  recorded  on  magnetic  tape.  The  resulting 
information  in  each  case  is  a continuous  time-history  record.  Prior  to 
each  test  shot,  a signal  that  corresponded  to  a known  value  of  the 
parameter  to  be  measured  was  recorded  on  the  parameter  data  channel. 

This  signal  is  hereafter  referred  to  as  the  "cal  step."  Parameters 
measured  during  a detonation  could  later  be  assigned  numerical  values  by 
comparison  with  the  cal  step. 

Filtering  and  digitizing  records 

U9.  Initial  reduction  of  the  time-history  records  from  magnetic 
tape  indicated  that  a substantial  amount  of  high-frequency  noise  was 
present  on  each  record.  The  magnitude  of  the  problem  was  considerably 
reduced  by  filtering  the  data  records  (including  cal  steps)  with  a 
10-kHz  low- pass  filter. 


50.  The  filtered  cal  steps  and  data  records  were  then  digitized 
at  1*8  kHz  and  rerecorded  on  magnetic  tape.  Using  an  editing  program 
(802-G9R0-101,  WES  computer  library),  the  desired  portions  of  the  cal 
step  and  corresponding  data  records  were  printed  from  the  magnetic  tape. 
The  printed  output  was  then  used  with  the  cal  step  value  to  determine 
the  conversion  factor  needed  for  each  tape  channel  to  change  the  digi- 
tized output  units  to  the  proper  units  of  measurement  recorded. 
Integration  of  data 

51.  The  filtered  digital  acceleration  and  velocity  data  were 
integrated  using  a previously  developed  program  amenable  to  digitized 
data  (803-G9RO-157,  WES  computer  library).  The  results  were  printed 

in  graphic  format;  sample  outputs  are  shown  in  Figures  22  and  23.  Fig- 
ure 22  shows  the  acceleration  data  from  an  accelerometer  mounted  in  the 
center  of  log  2,  shot  BB004;  the  integrated  velocity  and  displacement 
are  also  shown.  Figure  23  shows  the  results  of  measurements  from  a 
velocity  gage  in  the  center  of  log  2,  shot  BB003;  these  data  were  inte- 
grated to  give  displacement. 

Prediction  of  OP 

52.  The  initial  oscillograms  showing  OP  versus  time  were  diffi- 
cult to  interpret  and,  as  noted  above,  were  filtered  and  digitized  to 
aid  in  data  interpretation.  The  shape  of  the  resulting  traces  exhibited 
good  agreement  with  theory;  specifically,  an  almost  instantaneous  rise 
to  peak  OP  was  evident  followed  by  a rapid  initial  decay  in  the  general 
form  of  an  exponential  function  with  a negative  exponent.  The  final 
port  of  the  curve  was  characterized  by  a much  slower  (than  exponential) 
decay  until  zero  OP  was  attained. 

53.  Detailed  examination  of  the  initial  portion  of  the  traces, 
immediately  following  the  OP  rise,  showed  that  the  gages  tended  to 
overreact,  thereby  indicating  OP's  that  were  greater  than  the  ones  that 
actually  occurred.  This  problem  was  eliminated  by  following  the  gen- 
eral methods  developed  by  the  Ballistics  Research  Laboratory  for  esti- 

7 

mating  peak  OP.  This  technique  is  used  to  to  represent  the  OP  data  as 
a linear  function  of  time  and  to  estimate  peak  OP  at  the  beginning  of 
the  positive  phase  duration.  As  noted  in  paragraph  52,  the  final 
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Figure  22.  Velocity  and  displacement  obtained  from  digitized 
acceleration  data  (from  test  shot  BBOOU) 
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Figure  23.  Displacement  obtained  from  digitized  velocity  data 

(from  test  shot  BB003) 

portion  of  the  OP  decay  did  not  have  the  same  characteristic  envelope 
as  the  earlier  decay;  therefore,  the  last  two-thirds  of  the  decay  curve 
was  deleted  from  the  analysis.  The  remaining  data  were  represented 
analytically  in  the  following  manner. 

OP  = A * e"BT  (1) 

where 

OP  * overpressure,  psi 
A,B  * constants  for  given  data 
T * time,  msec 

»*7 


Taking  the  natural  logarithm  of  Equation  1 yields 


In  (OP)  « C - BT  (2) 

where  C » In  A . Equation  2 has  a linear  form.  At  T * 0 , the 
intercept  i3 


In  (OP)  = C 

If  an  instantaneous  rise  time  is  assumed,  the  constant  C and  the 
logarithm  of  the  peak  OP  are  the  same;  the  peak  OP  is  then 

OP  * eC 

This  method  of  estimating  the  peak  OP  represents  minimum  error  intro- 
duced by  the  initial  overreaction  of  the  gage. 

54.  Constants  B and  C were  determined  for  Equation  2 from  a 
least-squares  fit  to  the  data  points  in  the  initial  one-third  of  the 
pulse  duration.  Figures  24  and  25  show  the  raw  oscillogram  data,  digi- 
tized and  filtered  data,  and  the  OP  estimated  for  5.49  and  8.8  m, 
respectively,  from  GZ.  Comparison  of  the  estimated  OP  values  and  the 
predicted  values  is  presented  in  paragraphs  73-77. 

High-Speed  Photography 

55.  High-speed  motion  picture  records  were  made  of  each  shot. 
These  records  were  used,  when  possible,  to  determine  the  final  veloci- 
ties of  the  flying  logs  and  to  plot  the  visible  shock  front.  As  stated 
in  paragraph  11,  two  noninstrumented  logs  were  used.  The  ends  of  the 
logs  were  painted  in  alternating  black  and  white  quadrants.  This  tech- 
nique provided  a method  of  observing  rotation  around  the  longitudinal 
axis  of  the  log  and  the  displacement  of  the  center  of  the  log  during 
free  flight. 

56.  Log  motion  was  analyzed  with  a Vangard  Photographic  Data 
Reduction  Unit,  Model  M-l6c,  manufactured  by  Vangard  Instrument  Corp. 
Log  displacement  was  determined  by  recording  the  x-y  coordinates  of 
the  center  of  the  flying  log  relative  to  a known  fixed  point.  Timing 
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a.  OSCILLOGRAM  Of  STATIC  OVERPRESSURE  TRACE 


b.  OVERPRESSURE  TRACE  FILTERED  AND  DIGITIZED 


C.  NATURAL  LOGARITHM  Of  FILTERED  AND  DIGITIZED 
STATIC  OVER  PRESSURE  DATA 

Figure  25.  Estimation  of  peak  OP  at  8.8  m 
from  GZ 
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marks  on  the  film  in  conjunction  with  the  measured  coordinates  provided 
a time-history  of  the  log's  trajectory. 

57.  As  stated  in  paragraph  37,  three  plywood  backboards 
(1.2  by  2.k  m)  were  erected  at  the  test  site  to  study  the  shape  of  the 
shock  front  (Figure  13).  These  boards  were  positioned  perpendicular  to 
the  camera  and  were  painted  with  black  and  white  diagonal  stripes.  This 
configuration  and  the  variation  in  contrast  of  the  forest  in  the  back- 
ground provided  sufficient  image  distortion  to  locate  the  position  of 
the  shock  front  for  test  shots  BB002,  BB003,  BB006,  and  BB007.*  The 
shape  of  the  shock  front  was  analyzed  with  a Model  22U-A  Photo-Optical 
Data  Analyzer  manufactured  by  L-W  Photo,  Inc.  With  this  unit  individual 
frames  can  be  viewed  (Figure  26).  The  shape  and  position  of  the  shock 


Figure  26.  Sample  frame  of  high-speed  photographic 

record 


front  for  a given  frame  could  be  plotted  relative  to  GZ  by  rapidly 
viewing  adjacent  frames.  By  selecting  several  frames  during  the  time 
that  the  shock  front  was  in  the  field  of  view,  a succession  of  shock- 
front  shapes  could  be  plotted. 

* The  high-speed  photography  from  test  shots  BB001,  BBOOU,  and  BB005 
was  not  of  sufficient  quality  to  be  amenable  to  analysis. 
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PART  IV:  DATA  ANALYSIS  AND  RESULTS 


Quality  of  Test  Data 


58.  Prior  to  validation  of  the  theoretical  loading  function"*-  by 
using  results  from  the  test  shots,  it  was  necessary  to  examine  the  qual- 
ity of  the  data.  It  was  decided  that  the  following  criteria  must  be 
satisfied  before  the  data  could  be  considered  valid: 

a.  The  frame  of  reference  for  each  gage  measuring  displace- 
ment or  each  gage  dependent  on  displacement  for  genera- 
tion of  data  must  remain  fixed  with  respect  to  GZ  while 
being  acted  upon  by  the  shock  front  (this  included  all 
gages  except  the  velocity  gages  and  accelerometers 
mounted  on  the  logs).  This  frame  of  reference  could  be 
distorted  by  both  ground  and  air  shock. 

b.  The  shock  front  must  be  horizontally  symmetrical,*  i.e., 
there  must  be  no  Jetting  effects  (see  paragraphs  6U-66 
for  explanation  of  Jetting)  that  would  generate  abnor- 
mally high  OP's  for  a given  distance  from  GZ. 

c_.  The  data  must  be  taken  in  the  Mach  stem**  region.  This 
is  necessary  because  available  data  amenable  to  compari- 
son are  valid  only  in  this  region  of  the  shock  front. 

Also  the  loading  function  can  accommodate  cnly  one  rise 
and  decay,  i.e.,  above  the  Mach  stem  both  the  incident 
and  reflected  waves  pass  a given  point. 

d,.  The  yield  data  must  be  compatible  with  those  from  ac- 
cepted sources  of  blast  data. 

Stability  of  the  frame  of  reference 

59-  Moving  and  fixed  frames  of  reference  (with  respect  to  GZ) 
were  used  in  the  test  program.  The  flying  logs  instrumented  with  ve- 
locity and  acceleration  gages  used  a moving  frame  of  reference.  Exper- 
imental equipment,  the  measurement  capability  of  which  depended  on  a 
fixed  frame  of  reference,  included  the  displacement-time  apparatus, 

TOA  gages,  and  OP  gages.  As  mentioned  in  Part  II,  each  gage  was 
mounted  so  that  effects  caused  by  its  translation  or  vibration  with 
respect  to  its  frame  of  reference  were  minimized.  A gage  and  its 


* Circular  in  shape  as  viewed  in  the  horizontal  plane. 

**  The  shock  front  formed  by  the  fusion  of  the  incident  and  reflected 
shook  fronts  from  an  explosion;  see  Reference  6. 
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respective  frame  of  reference  must  remain  fixed  with  respect  to  each 
other;  otherwise,  measurements  will  be  with  respect  to  the  frame  of 
reference  of  the  gage  and  will  have  no  meaning  with  respect  to  the  fixed 
frame  of  reference  of  the  experiment,  i.e.,  with  respect  to  GZ. 

60.  Supporting  hardware  for  the  fixed  frame  of  reference  experi- 
mental apparatus  was  fabricated  according  to  one  of  two  techniques. 

a.  The  TOA  gages  (Figure  7)  and  some  of  the  OP  gages 
(Figure  6)  were  mounted  between  the  source  of  the  air- 
blast  and  the  support  hardware  so  that  the  gage  was 
acted  upon  by  the  shock  front  before  the  support  hard- 
ware. This  method  allowed  the  data  to  be  taken  before 
any  vibration  or  translation  of  the  support  hardware 
could  take  place. 

b.  The  superstructure  of  the  support  apparatus  for  each 
displacement-time  gage  (Figures  Ua  and  13)  was  designed 
to  withstand  any  effects  caused  by  the  incident  shock 
front.  Detailed  surveys  after  each  shot  showed  no  dis- 
placement of  any  rigid  support  apparatus  at  8.8  m from 
GZ.  Survey  accuracy  was  better  than  1.5  mm.  Tests 
BB006  and  BB007  were  arranged  so  that  the  distance  from 
GZ  was  less  than  8.8  m for  the  galvanized-pipe  support 
apparatus  (paragraphs  39  and  Uo).  No  displacement  or 
structural  damage  was  noted  after  shot  BB006  (OP  of 
200  psi);  however,  the  support  apparatus  did  fail 
during  shot  BB007  (OP  of  700  psi;  see  Figure  17). 

61.  Although  the  methods  of  fabricating  the  support  hardware  de- 
scribed above  were  effective  against  airblasts,  the  possibility  had  to 
be  considered  that  the  compression  wave  induced  in  the  earth  by  the  ex- 
plosion might  travel  ahead  of  the  shock  front  and  invalidate  mounting 
method  a,  paragraph  60,  i.e.,  vibration  of  the  aluminum  poles  used  to 
support  the  TOA  gages  might  break  the  pencil  lead  (paragraphs  18-21) 
before  the  shock  front  arrived,  or  vibration  of  the  stands  used  to  sup- 
port the  OP  gages  might  generate  erroneous  OP  data  (see  paragraph  17). 

To  investigate  the  magnitude  of  this  problem,  a seismic  refraction  sur- 
vey (paragraph  5)  was  conducted  prior  to  test  shot  BB001.  Results 
showed  that  the  compression  wave  velocity  in  the  near-surface  layer  was 
300  m/sec.  At  8.8  m from  GZ  (location  of  gages  for  most  tests),  the 
shock-front  velocity  for  a 0.09-metric-ton  yield  had  decayed  to  a 
predicted  value  of  6^3  m/sec;  therefore,  prior  to  initial  testing  it 
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seemed  probable  that  earth  motion  was  of  no  consequence.  This  assump- 
tion was  verified  by  the  earth-movement  canister  placed  at  8.8  m from 
GZ  for  tests  BB001  and  BB002  (see  paragraphs  23  and  3b).  For  test  BB001, 
the  TOA  gages  indicated  that  the  shock  front  arrived  at  8.8  m from  GZ 
at  6.123  msec  (average)  after  detonation;  the  earth-movement  canisters 
were  activated  no  earlier  than  8.0  msec  and  as  late  as  11.0  msec.  For 
the  second  test,  the  TOA  gages  indicated  that  the  shock  front  arrived 
at  8.8  m at  6.7  msec  (average)  after  detonation;  the  earth-movement 
canisters  were  activated  no  earlier  than  7.8  msec  and  as  late  as 
10.8  msec.  Therefore,  earth  motion  had  no  effect  on  the  gages  during 
the  period  when  data  were  being  taken. 

62.  Further  experiments  to  study  the  effects  of  earth  motion 
were  conducted  during  shot  BB005.  Three  TOA  gages  were  placed  in  cylin- 
drical cavities  as  described  in  paragraph  38;  following  the  detonation, 
the  gages  were  examined.  None  of  the  pencil  leads  in  the  TOA  gages 
were  broken,  indicating  that  earth  motion  did  not  affect  the  perfor- 
mance of  the  gages. 

63.  Based  on  the  results  of  the  investigation  discussed  above, 
it  was  believed  safe  to  assume  that  earth  motion  induced  by  the  explo- 
sions had  no  effect  on  the  quality  of  the  data  for  this  test  series. 
Shock-front 

symmetry  in  horizontal  plane 

64.  The  shock  front  resulting  from  a detonation  must  be  proven 
to  be  symmetrical  in  the  horizontal  plane  (circular)  so  that  all  data 
taken  at  the  same  distance  from  GZ  can  be  interpreted  with  equal  weight. 
The  major  phenomenon  that  destroys  the  symmetry  of  a shock  front  is  re- 
ferred to  as  "Jetting.'1 

65.  Generally,  Jets  occur  close  to  GZ  and  are  a result  of  some 
of  the  explosive  material  being  displaced  before  reacting.  At  some 
time  later,  the  detonation  catches  up  with  the  explosive  material  and 
produces  a Jet  (sometimes  luminous).  Associated  with  the  Jet  is  an 
OP  that  is  greater  than  would  normally  be  anticipated  at  that  distance 
from  GZ. 
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66.  The  effects  due  to  jetting  were  analyzed  by  evaluating 
arrival  times  at  different  azimuths  but  at  the  same  distance  from  GZ. 
Equal  arrival  times  indicated  that  the  shock  front  was  symmetrical  in 
the  horizontal  plane,  and  there  vas  no  detectable  Jetting.  Arrival- 
time data  were  available  from  two  sources,  the  TOA  gages  and  the  acti- 
vation time  of  the  OP  gages.  Tables  6 and  7 show  a complete  compila- 
tion of  TOA  and  OP  gage  data,  respectively,  for  the  entire  test  series. 
All  TOA  data  in  which  there  were  at  least  two  measurements  at  the  same 
radial  distance  from  GZ  were  complied  from  Table  6 and  listed  in 
Table  6.  A similar  compilation  was  not  prepared  from  Table  7 because 
of  the  variable  response  times  of  the  OP  gages;  however.  Table  7 is 
presented  for  completeness  and  general  interest. 

67.  Table  8 contains  the  TOA  data  for  test  shots  BB001  and 
BB002.  For  test  shot  BB001,  four  TOA  gages  were  spaced  90  deg  apart, 
all  at  8.8  m from  GZ.  Arrival  times  of  6.1,  6.2,  6.2,  and  6.0  msec 
were  recorded.  The  small  differences  in  arrival  time  indicate  that 
symmetry  was  almost  achieved  in  the  horizontal  plane  at  8.8  m.  These 
differences  should  be  considered  within  the  limits  of  experimental 
error. 

68.  Symmetry  in  the  horizontal  plane  was  studied  in  greater 
detail  during  the  second  test,  BB002  (see  paragraph  35).  In  addition 
to  the  four  TOA  gages  used  in  BB001,  two  strings  of  gages  were  added 
to  the  test  configuration.  These  strings  were  perpendicular  to  each 
other.  The  TOA  gages  were  spaced  such  that  the  estimated  travel  time 
of  the  shock  front  from  one  gage  to  the  next  would  be  2 msec.  The 
results  are  presented  in  Table  8.  The  symmetry  distortion  was  again 
evaluated  by  comparing  arrival  times  for  gages  at  the  same  radial 
distance  from  GZ.  A 10  percent  variation  in  TOA  is  not  unreasonable 
to  expect  for  two  gages  placed  at  the  same  distance  from  GZ  for  a 
given  shot.  The  data  in  Table  8 meet  this  criterion,  indicating  that 
the  shock  front  propagated  away  from  GZ  in  a nearly  symmetrical 
manner . 
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Shape  of  the  shock  front 
in  the  vertical  plane  and 
the  path  of  the  Mach  stem 

69.  Several  experiments  were  included  in  the  test  program  to 
study  the  shape  of  the  shock  front  in  a vertical  plane  passing 
through  GZ  and  the  path  of  the  Mach  stem.  Data  on  the  shape  of  the 
shock  front  in  a vertical  plane  were  required  if  it  became  apparent 
that  the  plane-wave  shock  front  used  in  the  original  loading  function 
did  not  correctly  simulate  the  actual  shock  front.  Study  of  the  Mach 
stem  path  was  necessary  to  verify  that  the  data  vere  taken  in  the  Mach 
stem  region;  data  taken  above  the  Mach  stem  could  not  be  scaled  to 
other  yields,^  i.e.,  the  scaling  laws  are  valid  only  in  the  Mach  stem 
region.  In  addition,  the  loading  function  is  valid  only  in  the  Mach 
stem  region. 

70.  A shock  wave  generated  by  an  explosion  in  air  is  approxi- 
mately spherical.  If  the  explosion  is  close  to  the  ground,  a Mach  stem 
develops  as  soon  as  the  shock  wave  strikes  the  ground.  Since  the  veloc- 
ity of  the  Mach  stem  away  from  GZ  is  slightly  greater  than  that  of  the 
incident  shock  wave,  the  shape  of  the  shock  wave-Mach  stem  combination, 
as  viewed  in  a vertical  plane  passing  through  GZ,  is  a complex  curve, 
with  a breakpoint  (called  the  triple  point)  at  the  junction  between  the 
Mach  stem  and  the  incident  shock  wave  (Figures  27  and  28). 

71.  Experiments  to  examine  the  inclination  of  the  shock  front 
were  conducted  during  shots  BB001  and  BB002  (see  paragraphs  33  and  3M. 
Pairs  of  TOA  gages  were  placed  at  the  same  distance  from  GZ,  but  at  dif- 
ferent elevations  (see  Table  U and  Figure  7b).  For  the  shock  front  to 
be  spherical  in  nature,  the  lower  TOA  gage  of  a pair  should  be  acti- 
vated before  the  upper  gage.  The  resulting  data  are  presented  in 
Table  9.  The  majority  of  the  data  were  lost  in  shot  BB002  because  of 
circuit  failure.  The  three  pairs  of  data  that  are  valid  for  each  shot 
show  that  the  lower  TOA  gages  were  indeed  activated  first,  indicating 
that  the  shock  front  was  inclined  toward  GZ. 

72.  Photography  suitable  for  analysis  of  the  shape  of  the  shock 
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«.  SHOT  NO.  BB002  - OAT  A OBTAINEO  FROM  BLACK  AND 
WHITE  FILM  EXPOSED  AT  7BSB  FRAMES  PER  SECOND 


NOT  E 

1 ALL  HORIZONTAL  DISTANCES  FROM  CZ  ARE 
ACCURATE  IN  RLAN  VIEW,  RUT  ARE  ONLY 
APPROXIMATE  IN  THE  FIELDS  OF  VIEW  SHOWN 

2 FRAME  NO'S  REFER  TO  THE  NO  OF  « MM  FILM 
FRAMES  EXPOSED  SINCE  DETONATION. 


b.  SHOT  NO.  BB003  ■ DATA  OBTAINED  FROM  COLOR  FILM 
EXPOSED  AT  2000  FRAMES  PER  SECOND 

Figure  27.  General  shape  of  shock  front  for  test  shots 

BB002  and  BB003 
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«.  SHOT  NO.  BBOOfi  - DATA  OBTAINED  FROM  COLOR  FILM 
EXPOSED  AT  S2I0  FRAMES  PER  SECOND 


NOTE 

I ALL  HORIZONTAL  DISTANCES  FROM  CZ  ARE 
ACCURATE  IN  FLAN  VIEW,  BUT  ARE  ONLY 
APPROXIMATE  IN  THE  FIELDS  OF  VIEW  SHOWN. 

Z.  FRAME  NO'S  REFER  TO  THE  NO  OF  It  MM  FILM 
FRAMES  EXPOSED  SINCE  DETONATION 


b.  SHOT  NO.  BB007  - DATA  OP  SAINED  FROM  COLOR  FILM 
EXPOSED  AT  SISO  FRAMES  PER  SECOND 


Figure  28.  General  shape  of  shock  front  for  test  shots 

BB006  and  BB007 
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front  was  obtained  for  test  shots  BB002,  BB003,  BB006,  and  BB007.*  The 
framing  rates  were  7859,  2000,  5210,  and  5150  f/s,  respectively.  The 
shock- front  shapes  were  constructed  using  the  methods  described  in  para- 
graph 57.  The  shapes  of  the  shock  fronts  as  they  propagated  away  from 
GZ  for  shots  BB002,  BB003,  BB006,  and  BB007  are  shown  in  Figures  27  and 
28.  The  path  of  the  triple  point  is  sufficiently  obvious,  at  least  in 
Figure  28,  to  conclude  that  all  measurements  of  shock-front  T0A,  OP, 
ground  motion,  and  log  motion  were  obtained  in  the  Mach  stem  region. 

Compatibility  of  yield 

data  with  accepted  source  data 

73.  The  experimental  yield  data  were  compared  with  the  following 
accepted  data  sources: 

ti.  Ballistics  Research  Laboratory  (BRL)  OP  tost  data,  which 
were  scaled  from  3.6-  to  0.09-metric-ton  yield. 

]a.  Air  Force  Weapons  Laboratory  (AFWL)  test  data,**  which 
were  scaled  from  normalized  data. 

£.  Naval  Ordnance  Laboratory  (NOL)  test  data,'’  which  were 
scaled  from  normalized  data. 

71*.  These  sources  were  selected  for  comparison  for  the  following 
reasons: 

a.  The  BRL  data  were  scaled  down  from  data  that  were  in 
good  agreement  with  data  acquired  at  Fort  Benning.^ 

After  the  Fort  Benning  test,  the  BRL  data  were  used  to 
generate  the  generalized  stem-failure  criterion  and  to 
structure  the  loading  function. 

b.  The  AFWL  test  data  were  acquired  from  the  detonation  of 
a 0.U5-metric-ton  spherical  charge  of  TNT,  with  the 
height  of  burst  being  b.57  m;  U . 57  m scales  down  to 
1.33  m for  0.09-metric-ton  TNT  yield.  This  value 
(1.33  m)  is  very  close  to  the  height  of  burst  used  for 
the  test  series  described  herein  (0.72  m). 

£.  The  NOL  test  data  were  valid  for  a spherical  charge  in 
free  space,  i.e.  so  far  above  the  ground  that  no  Mach 
stem  forms.  These  data  were  used  only  to  compare  the 
shock-front  characteristics  of  a free-space  explosion 
with  those  of  an  explosion  at  a finite  height  above 
ground  where  a distinct  Mach  stem  was  formed. 


* The  high-speed  photography  from  test  shots  BB001,  BBOOU,  and  BB005 

was  not  of  sufficient  quality  to  be  amenable  to  analysis. 
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75.  Peak  OP  and  TOA  for  each  of  the  above  sources  were  scaled  to 
0.09-metric-ton  TNT  yield.  The  scaled  values  are  plotted  in  Figures  29 
and  30,  respectively.  The  WES  TOA  and  OP  data  are  in  good  agreement 
with  the  data  from  the  other  sources.  The  compatibility  of  data  on  both 
plots  indicates  that  the  desired  energy  yields  were  obtained  from  the 
AN/FO  explosives,  and  that  existing  cube-root  scaling  laws  may  be  ap- 
plied with  confidence  to  yields  as  low  as  0.09-metric-ton  TNT. 

76.  The  NOL  data  appear  to  produce  a lower  OP  and  later  arrival 
time  than  the  other  available  data.  This  is  because  the  other  data  are 
measured  in  the  Mach  stem  region,  and  the  NOL  data  are  valid  only  for 
free  space,  i.e.,  the  Mach  stem  is  the  composite  of  both  the  incident  and 
reflected  waves,  whereas  the  NOL  data  result  from  the  incident  shock 
front  only.  The  major  differences  appear  to  be  close  to  GZ;  however,  at 
large  distances  from  GZ,  the  differences  become  much  smaller.  The  BRL, 
AFWL,  and  NOL  data  shown  in  Figures  29  and  30  were  standardized  such 
that  they  are  independent  of  temperature  and  barometric  pressure. 

77.  Numerical  comparison  of  WES  and  BRL  TOA  and  OP  data  are  shown 
in  Tables  10  and  11,  respectively.  These  data  were  not  corrected  for 
variations  in  meteorological  conditions.  Variation  in  the  barometric 
pressure  is  insignificant  for  the  data  shown  in  both  tables.  Tempera- 
ture introduces  insignificant  error  at  8.8  m and  beyond.  The  available 
data  indicate  that  approximate  symmetry  of  the  shock  front  was  achieved 
in  the  WES  test  shots  at  about  5*5  m.  However,  the  fireball  normally 
extended  out  to  about  8.0  m (Figure  31). 

78.  Existing  methods  of  predicting  TOA  and  OP  are  valid  only  when 
the  available  energy  is  entirely  kinetic,  i.e.,  outside  the  fireball. 

This  kind  of  correction  cannot  be  applied  inside  the  fireball,  since  in 
that  situation  chemical  reactions  are  continuing  and  adding  heat  energy 
to  the  expanding  gases.  The  kinetic  instabilities  resulting  from  this 
situation  make  it  invalid  to  apply  normal  temperature  corrections  to  the 
region  inside  the  fireball,  which  in  the  case  of  the  WES  tests  implies 
any  distance  less  than  about  8.0  m.  Thus,  even  though  the  shock  front 
appears  to  have  attained  symmetry  at  about  5-5  m (which  is  generally 
regarded  as  the  distance  from  GZ  beyond  which  normal  TOA  and  OP 


60 


Figure  29.  Peak  OP  versus  distance  from  GZ  for  selected  data 


ARRIVAL  TIME 


oi stance  from  ground  zero.  n. 


Figure  30.  Arrival  time  versus  distance  from  GZ  for  selected  data 
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Figure  31.  Explosion  of  test  shot  BB002.  Fireball  has 
expanded  to  maximum  value  of  approximately  8.0  m 

prediction  methods  could  be  re  .iably  used),  the  TOA's  and  OP's  in  the 
annulus  between  5.5  and  8.0  m cannot  be  reliably  predicted. 

79-  Of  the  56  TOA  measurements  obtained  at  distances  from  8.8  m 
up  to  22.5  m,  89  percent  are  within  10  percent  of  the  values  determined 
from  the  BRL  data  (Table  10).  Of  the  17  peak  OP  values  determined  from 
the  WES  experiments,  82  percent  are  within  15  percent  of  peak  values 
obtained  from  BRL  data,  and  100  percent  are  within  21  percent  of  the 
BRL  values  (Table  11).  A 20  percent  variation  in  peak  OP  measurement 
is  not  uncommon  for  similar  gages  placed  at  the  same  distance  from  GZ 
for  a given  shot. 

80.  To  evaluate  the  compatibility  of  the  WES  TOA  data  and  OP 
data,  the  TOA  data  were  converted  to  OP  data  as  follows:  The  distance 

from  GZ  was  plotted  versus  TOA  using  all  valid  data.  A least-squares 
polynomial  equation  was  generated  using  these  data.  The  derivative  of 
this  equation  yields  the  instantaneous  velocity  of  the  shock  front  for 
a given  arrival  time.  This  arrival  time  corresponds  to  a given  dis- 
tance from  GZ;  hence,  the  shock- front  velocity  can  be  plotted  as  a func- 
tion of  distance  from  GZ.  The  Rankine-Hugoniot  equation  for  shock-wave 
velocity  (V)  in  the  Mach  stem  region  is^ 


where 


. ■ ,(■  • C 


(3) 


Cq  = velocity  of  sound  at  ambient  temperature 
P = peak  overpressure 
Pq  = ambient  atmospheric  pressure 

Rearranging  Equation  3 yields 


Therefore,  a curve  of  peak  OP  versus  distance  from  GZ  can  be  plotted 
(distance  from  GZ  is  a function  of  velocity  as  stated  above).  This 
curve  is  shown  in  Figure  29  along  with  measured  WES  data,  and  BRL,  AFWL, 
and  NOL  peak  OP  data.  The  curve  appears  to  be  in  good  agreement  with 
the  other  curves. 

Summary 

81.  The  data  obtained  during  the  WES  test  series  agree  with 
theory  within  anticipated  instrumental  and  observational  error  based  on 
the  conditions  specified  in  paragraph  58,  i.e.,  the  frame  of  reference 
of  each  sensor  was  not  influenced  by  the  incident  shock,  the  shock  front 
exhibited  symmetrical  propagation  characteristics,  all  data  were  taken 
in  the  Mach  stem  region,  and  yield  data  from  the  WES  test  series  are 
compatible  with  yield  data  from  accepted  sources  of  blast  data. 

Validation  of  Tree  Stem  Loading  Function 

82.  The  major  objective  of  the  test  program  was  to  validate  the 
theoretical  loading  function  used  in  the  solution  of  the  partial  differ- 
ential equation  (PDE)  describing  the  dynamic  deformation  of  a right- 
circular  cylindrical  beam,  i.e.  the  dynamic  motion  of  an  individual 
tree  stem  being  acted  upon  by  a shock  front. ^ The  solution  of  the  PDE 
was  required  in  order  to  relate  the  geometric  stem  curvature  required 
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for  stem  failure  to  the  height  on  the  stem  (remnant  height)  at  which 
failure  occurs.  Development  of  the  stem  failure  criterion  and  solution 
of  the  PDE  are  discussed  in  Appendix  A.  Using  the  failure  criterion 
and  the  PDE  solution,  a nomographic  technique  was  developed  for  estimat- 
ing stem-remnant  height  at  a given  distance  from  GZ  for  wide  variations 
in  stem  properties  and  weapon  yields.  The  development  of  this  procedure 
is  described  in  Appendix  B. 

83.  The  tree  stem  discussed  above  is  represented  in  this  valida- 
tion by  the  logs  used  in  this  test  program  (see  paragraph  9).  The 
loading  function  was  validated  by  considering  the  momentum  transfer  per 
unit  length  of  the  log  between  the  shock  front  and  the  log  during  the 
positive  pulse  duration  of  the  shock  front,  i.e. 


M = 


F(x,t)  dt 


(h) 


where 

M = momentum  transfer  per  unit  length  of  the  log 

Tp  = length  of  positive  pulse  duration 

F(x,t)  = loading  function  that  characterizes  the  dynamic  loading 

of  the  log  by  the  shock  front  (see  Appendix  A for  deriva- 
tion of  loading  function) 

Integration  of  Equation  U by  expressing  F(x,t)  as  a change  of  momen- 
tum with  respect  to  ti^ie  yields 


(5) 


where 

m = mass  of  log 
£ = length  of  log 

Vj,  = velocity  of  log  at  end  of  positive  pulse  duration  (referred 
to  herein  as  final  velocity) 

= velocity  of  log  at  beginning  of  positive  pulse  duration 
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was  always  zero  for  the  WES  test  series  because  the  logs  were  at 
rest  when  the  shock  front  initially  came  in  contact  with  them.  Experi- 
mental data  for  V^,  are  available  from  the  field  procedures  previously 
discussed  in  paragraphs  1+7,  51,  55,  and  56;  also  see  paragraphs  86-91. 

A theoretical  prediction  for  V can  be  obtained  by  integrating  Equa- 
tion 4,  using  the  derived  loading  function;  the  resulting  value  of  M 
can  be  equated  with  Equation  5 to  obtain  the  desired  theoretical  value 
of  Vf  . Agreement  of  the  experimental  and  theoretical  values  of 
for  the  same  data  will  validate  the  loading  function. 

Theoretical  determination  of  V 

84.  Integrating  Equation  4 using  the  derived  loading  function 
(see  Appendix  A)  shows  that  momentum  transfer  per  unit  length  is 

M = 2rJT  (6) 

where 

r = radius  of  the  stem 

JT  = the  total  effective  pressure  impulse;  see  paragraphs  29 

through  38  of  Reference  1 for  discussion  of  the  total  effec- 
tive pressure  impulse 

J^,  can  be  calculated  from 

JT  * JCd  * Pr(r) 

where 

J » dynamic  pressure  impulse 
* drag  coefficient  of  stem 
Pr  * reflected  pressure 
U * velocity  of  shock  front 
JC^  is  generated  by  the  air  drag  on  the  stem  within  the  shock  front; 

Pr  is  due  to  a pressure  differential  generated  between  the  side  facing 
the  incident  shock  front  and  the  lee  side  of  the  tree  stem;  and  2r/U 
is  the  shock-front  transit  time  across  the  stem.  The  derivation  of 


! 

(7) 
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Equation  7 is  shown  in  Reference  1 on  pages  25  and  26.  The  last  term 
in  Equation  7 can  be  expressed  in  terms  of  known  values  of  the  dynamic 
pressure  impulse  (J0)  as  follows 

Pr(r)  - °'WB  (8> 

For  a stem  radius  of  10  cm,*  the  term  Pr(2r/U)  was  calculated  using 
OP  values  from  Figure  8 of  Reference  1.  The  results  indicate  that  for 
a given  distance  from  GZ,  the  calculated  value  of  Pr(2r/U)  is  equal 
to  approximately  one-half  the  corresponding  dynamic  pressure  impulse 
(also  shown  in  Figure  8).  To  satisfy  any  stem  diameter,  both  terms  in 
Equation  8 can  be  multiplied  by  the  term  r/10  . According  to  cube- 
root  scaling  laws,  the  dynamic  pressure  impulse  can  be  computed  for  any 
weapon  yield  (W)  for  a given  distance  from  GZ  if  the  value  for  the 
dynamic  pressure  impulse  ( J_ ) is  known  for  another  yield  (W_ ) at  the 
same  distance  from  GZ,  i.e. 


J = 


(9) 


Substituting  Equations  8 and  9 into  Equation  7 and  the  resulting  value 
of  J^,  into  Equation  6 yields 


M = d 


JB  X 


(10) 


where  d = 2r  . Equating  Equations  5 and  10  and  solving  for 
yields  the  desired  theoretical  final  velocity: 


V 


f 


di 

m 


(11) 


* Arbitrarily  selected  value;  any  real  number  can  be  used. 
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85 . A sample  calculation  to  obtain  a theoretical  value  of 
(Equation  11)  is  shown  below  using  typical  physical  properties  of  the 
logs  used  in  the  WES  test  series  (compiled  in  Table  12).  Consider  shot 
BB001,  log  1: 

m « 9054  g 
l » 111.8  cm 
d ■ 10  cm 

Also 

2 

Jg  * 2434  dyne  * sec/cm  (predicted  dynamic  pressure  impulse  at 
8.8  m from  GZ  for  a 0. 09-metric-ton  yield) 

For  all  calculations  of  V^, 

W « 0.09  metric  ton 
Wg  * 3.6  metric  tons* 

Cd  * 0.75** 

Then 


= 481.8  cm/sec 


Experimental  determination  of  V^, 

86.  The  final  velocity  (V^.)  was  obtained  experimentally  from  13 
logs.  As  stated  in  paragraph  12,  each  of  the  logs  was  111.8  cm  long, 
was  either  10  or  20  cm  in  diameter,  and  varied  in  mass  from  9*054  to 
36,083  g.  Twelve  had  rough,  natural  surfaces  and  one  had  a smooth  sur- 
face. The  logs  were  instrumented  to  measure  displacement,  velocity,  or 
acceleration.  Some  logs  were  instrumented  to  measure  more  than  one  of 


* See  Figure  8 of  Reference  1. 

**  Predicted  drag  coefficient  for  a cylinder  with  a surface  similar  to 
that  of  a log. 
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the  parameters.  The  displacement  measurements  vere  taken  from  three 
different  sources:  (a)  shot  BB001,  log  1 and  log  2 data  were  obtained 

from  the  displacements  measured  on  high-speed  movies  of  the  blast;  (b) 
shot  BB002 , log  1 data  were  taken  from  the  wire-breaker  device;  and  (c) 
the  displacement  data  for  shots  BB003-BB006  were  taken  from  the  magnetic 
wire-displacement  apparatus.  The  derivative  of  the  displ  icement-time 
data  yielded  the  desired  final  velocity.  Velocity  data  were  also  ob- 
tained from  direct  output  of  the  velocity  gages.  The  acceleration 
records  were  integrated  when  possible  to  obtain  velocity  data. 

87.  Eleven  of  the  logs  from  which  experimental  velocity  data  were 
obtained  were  statically  positioned  at  8.8  n from  GZ  prior  to  detona- 
tion; one  was  at  12.2  m,  and  one  was  at  4.1*  m (paragraphs  33-40).  These 
data  for  each  log,  plus  the  theoretically  calculated  and  experimentally 
measured  horizontal  velocities,  are  shown  in  Table  12.  Experimental 
data  are  shown  for  gages  positioned  in  the  center  and  in  both  ends  of 
the  log.  A value  that  can  be  interpreted  as  the  representative  hori- 
zontal velocity  of  the  center  of  mass  of  the  log  was  calculated  by  using 
simple  or  weighted  averages.  The  averages  were  obtained  by  classifying 
the  horizontal  velocity  data  into  one  of  the  following  categories: 

a.  One  or  more  values  at  only  one  of  the  three  gage 
positions. 

b.  One  or  more  values  at  both  end  mounting  positions,  but 
not  at  the  center. 

£.  One  or  more  values  at  each  of  the  three  mounting 
positions. 

The  representative  horizontal  velocity  of  the  center  of  mass  of  the  log 
was  then  obtained  in  the  following  manner:  Category  a data  were  aver- 

aged. A value  was  obtained  for  category  b data  by  averaging  all  data 
at  each  end,  and  then  averaging  the  resulting  end  values.  A weighted 
value  was  obtained  from  category  £ data  by  averaging  the  data  at  each 
position.  The  end  values  were  averaged;  the  result  of  this  average 
was  averaged  with  the  average  center  value  to  obtain  the  representative 
velocity. 

88.  Initial  inspection  of  the  measured  horizontal  velocities  for 
the  12  logs  with  natural  surfaces  (Table  12)  indicates  the  following: 
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a.  The  neasured  horizontal  velocity  was  somewhat  higher 
when  obtained  from  a velocity  gag''  than  from  an  acceler- 
ometer gage  when  data  are  compared  from  gages  mounted  at 
the  same  position  o\  the  log. 

b.  The  horizontal  velocity  appears  to  be  different  at  oppo- 
site ends  of  logs  that  had  instrumentation  on  both  ends. 

£.  The  difference  in  the  log  velocities  at  either  end  of  the 
log  appears  to  be  minimized  when  displacement-time  gages 
are  used. 

89.  The  data  in  Table  12  indicate  seven  instances  in  which  veloc- 
ity gages  measured  higher  values  than  those  measured  by  accelerometers 
for  the  final  log  velocity.  Even  though  this  is  not  enough  data  to 
rigorously  analyze  the  problem,  some  speculation  is  possible.  Errors 
due  to  integration  or  an  insufficient  sample  of  the  acceleration-time 
record  were  initially  eliminated.  Errors  due  to  nonrigid  mounting  were 
possible,  i.e.,  the  gages’  frames  of  reference  were  not  fixed  with  re- 
spect to  one  another;  however,  the  most  probable  cause  of  the  higher 
values  is  the  rotation  of  the  log  disturbing  the  calibration  of  the 
velocity  gages. 

90.  Rotation  of  the  log  was  noted  (high-speed  photography)  about 
the  vertical  and  lateral  axer.  of  the  log,  but  not  about  its  longitudinal 
axis.  Rotation  about  the  two  axes  was  probably  due  to  nonsymmetrical 
loading  of  the  log  (variation  in  OP  across  the  face  of  the  log),  varia- 
tion in  the  drag  coefficient  on  the  loading  surface,  or  the  center  of 
mass  of  the  log  being  different  from  the  geometric  center  of  the  log. 

The  rotational  motion  induced  by  one  or  more  of  the  above  parameters 
most  likely  accounts  for  the  differences  in  the  horizontal  velocities. 

91.  The  displacement-time  gages  produced  the  most  consistent  re- 
sults (see  Table  13)  i.e.,  the  horizontal  end  velocities  agree  more 
closely  than  those  velocities  determined  from  accelerometer  and  velocity 
data.  The  possibility  should  be  considered  that  the  drag  of  the  magne- 
tized wire  on  the  Teflon  guide  (paragraph  28)  reduced  the  difference  in 
the  end  velocities.  The  coefficient  of  kinetic  friction  for  Teflon  is 
close  to  zero,  thus  drag  of  the  magnetized  wire  in  the  Teflon  was  prob- 
ably insignificant  and  was  assumed  to  be  so  in  computing  final  log 
velocities.  The  quality  of  data  generated,  plus  low  cost  and  ease 
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of  installation  (as  compared  with  velocity  gages  and  accelerometers 
mounted  in  the  logs),  indicates  that  a displacement-time  apparatus  is 
the  optimum  configuration  for  obtaining  impulse  data  from  free-flying 
logs. 

Comparison  of  theo- 
retical  and  experimental  data 

92.  Figure  32  is  a comparison  of  theoretically  predicted  and 
experimentally  measured  final  velocities  for  13  logs  put  into  free 
flight  by  impulse  loading  from  a passing  shock  front.*  The  predicted 
values  were  computed  with  Equation  11.  The  linear  relation  shown  on  the 
graph  indicates  that  the  measured  velocities  are  somewhat  higher  than 
the  predicted  velocities  (note  that  the  slope  of  the  line  is  not  unity). 
However,  if  a higher  value  were  used  for  the  drag  coefficient,  the  pre- 
dicted values  could  be  adjusted  for  a match.  The  predicted  velocity  for 
the  smooth  log  is  too  high.  Therefore,  the  drag  coefficient  would  have 
to  be  lowered  to  adjust  for  a match.  Thus,  the  drag  coefficient  used 
to  calculate  the  predicted  final  velocities  (0.75)  shown  in  Figure  32 
is  too  low  for  logs  with  a rough  surface  and  too  high  for  a log  with 
a smooth  surface.  If  it  is  assumed  that  the  pine  logs  are  rough  compared 
to  the  majority  of  tree  species,  0.75  may  be  a good  compromise.  The 
large  cluster  of  points  in  Figure  32  represents  data  collected  on  logs 
of  two  diameters  when  the  incident  OP  was  ho  psi.  Also  single  tests 
were  conducted  at  20  and  200  psi  as  indicated  on  Figure  32.  The  re- 
sulting data  appear  to  be  in  good  agreement  with  the  linear  fit,  i.e., 
the  20-psi  data  point  falls  directly  upon  the  linear  fit,  and  the  200- 
psi  data  point  slightly  above.  Although  there  is  some  discrepancy 
between  the  measured  and  predicted  values  of  V^,  , the  data  spread  in 
Figure  32  is  small  enough  to  conclude  that  the  stem  loading  function 
(paragraph  8)  used  in  the  tree  height  prediction  model  properly  repre- 
sents the  physical  phenomena. 


* Data  from  last  two  columns  of  Table  12. 
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AVERAGE  MEASURED  FINAL  VELOCITY,  CM/SEC 


Figure  32.  Comparison  of  predicted  and  measured  final 
velocities  for  logs 


PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


93.  The  following  conclusions  are  based  on  the  measurement  tech- 
niques developed  for  the  test  program  and  the  resulting  data: 

a.  Reliable  equipment  and  data  reduction  techniques  were 
developed  to  evaluate  the  impulse  loading  experienced 
by  a tree  stem  (flying  log)  being  acted  upon  by  a 
shock  front. 

b.  The  reduced  data  from  the  test  series  indicate  that  the 
loading  function  used  in  the  partial  differential  equa- 
tion describing  tree  stem  motion  under  the  influence  of 
a shock  front*  is  valid  for  the  data  presented. 


Recommendations 


9^.  Based  on  the  findings  of  this  study , it  is  recommended  that: 

a.  Comprehensive  tests  be  initiated  to  acquire  additional 
impulse  data  for  tree  stem  loading  for  OP's  greater  than 
and  less  than  kO  psi. 

b.  Additional  tests  be  conducted  to  determine  the  relative 
influence  of  the  drag  coefficient  on  the  theoretical 
loading  function  as  compared  with  other  parameters. 
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Table  1 


Equipment  Used 

in  Test 

Shots 

Shot  BBOO- 

Equipment 

1 

2 

1 i 2 

6 

I 

Noninstrumented  logs 

X 

X 

Instrumented  logs 

(See  Table  2) 

Explosive  container 

c 

c 

css 

s 

s 

Recorders 

X 

X 

XXX 

X 

X 

Overpressure  (OP)  gages 

(See  Table  3) 

Time-of-arrival  (TOA)  gages 

(See  Table  U) 

Questar  telescope 

X 

X 

XXX 

X 

X 

Earth-movement  canister 

X 

X 

Accelerometers 

(See  Table  2) 

Velocity  gages 

(See  Table  2) 

Displacement-time  gages 

(See  Table  2) 

Cameras 

Milliken 

X 

X 

X 

X 

X 

Kodak 

X 

X 

Nova 

X 

X 

X 

X 

Fairchild 

XXX 

X 

X 

Fastex 

X X 

Image  distortion  backboards 

X X 

X 

X 

Note:  x * apparatus  used;  c * cubical;  and  s * spherical. 


Table  2 

Inatruaented  Logi  Used  In  Tc»t  Shots 
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Table  3 

Overpressure  Gages  Used  in  Test  Shots 


Height  of  Gage 

Overpressure 

Distance  from 

Above 

Test 

Gage 

Azimuth* 

Ground  Zero 

Reference  Plane** 

Shot 

Designation 

deg 

m 

m 

1 

1 

90 

8 

.8 

-0, 

.71* 

2 

180 

3 

270 

! 

it 

23 

1 

0.76 

2 

1 

93 

5.1*9 

-0.91* 

2 

183 

5.1*9 

I 

3 

273 

8.8 

1 

1* 

23 

8 

.8 

0.56 

3 

1 

93 

5 

.1*9 

-0.60 

2 

183 

-0.60 

3 

273 

0, 

.0 

U 

1 

138 

2 

183 

3 

273 

1 

5 

1 

108 

-0! 

60 

2 

293 

3 

207 

8 

8 

1 

it 

3th 

0:0 

5 

331 

6 

303 

12, 

.2 

6 

1 

108 

5.1*9 

2 

292 

5.1*9 

3 

207 

8, 

.8 

7 

1 

108 

5.65 

2 

292 

6. 

.19 

3 

207 

10.7 

1 

• Az 

inuth  relative  to 

point  in 

center  of  large  concrete  block  which 

was  base  of  wi re-breaker  displacenent-time  stand  used  in  test  shots 
PB001  and  BB002.  Origin  of  coordinate  system  at  ground  zero. 

**  The  reference  plane  is  the  horizontal  plane  which  includes  the  point 
source  of  the  explosion  (72  cm  above  original  ground  surface).  Posi- 
tive values  indicate  height  above  the  plane  and  negative  below  the 
plane.  This  plane  was  the  same  for  each  test. 


Table  1* ** 

Time-of-Arrlval  Gages  Used  In  Test  Shots 


Height  of  Height  of 

Distance  Gage  | Distance  Gage 


* Azimuth  relative  to  point  in  center  of  large  concrete  block  which  was  base  of  wire-breaker 
displacement-tine  stand  used  in  test  shots  BB001  and  BB002.  Origin  of  coordinate  system  at 
ground  zero. 

**  The  reference  plane  i3  the  horizontal  plane  which  includes  the  point  source  of  the  explosion 
(72  cm  above  original  ground  Burface).  Positive  values  indicate  height  above  the  plane  and 
negative  below  the  plane.  This  plane  was  the  sane  for  each  test. 


Table  5 

Meteorological  Data  For  Shots  BB001  through  BB007 


Test 

Temperature 

Barometric 

Shot 

°C 

Pressure,  psi 

BB001 

16.7 

li».8 

BB002 

29.  k 

1U.5 

BB003 

32.8 

1U.7 

BBOOU 

30.0 

lU.7 

BB005 

3k. 7 

Ik. 6 

BB006 

30.6 

Ik. 7 

BB007 

35. U 

lk.7 

< 


Table  6 

Data  from  Tlmc-of-Arrlval  (TOA)  Gages 


Height  of  Gage 

Distance  from 

Above 

Arrival 

Test 

TOA  Gage 

Azimuth 

Ground  Zero 

Reference  Plane** 

Time 

Shot 

Designation 

deg* 

m 

m 

msec 

BB001 

9 

0 

8.8 

0 

0 

6.1 

10 

90 

8.8 

6.2 

1 

180 

8.8 

6.2 

8 

270 

8.8 

6.0 

2 

180 

10.8 

8.6 

5 

10.9 

O 

37 

10.2 

3 

12.3 

0 

0 

12.1 

6 

12.3 

2 

37 

ll*.0 

1* 

13.8 

0. 

0 

15.0 

7 

13.8 

o 

37 

16.8 

BB002 

53 

0 

8.8 

0. 

0 

6.7 

51 

270 

8.8 

0. 

0 

6.7 

2U 

180 

10.1*5 

1. 

00 

ll.l* 

25 

180 

11.6l 

1. 

00 

13.1 

26 

180 

12.65 

1. 

00 

15.2 

27 

180 

13.65 

1 

00 

17.1* 

itl 

90 

1U.56 

0. 

0 

17.2 

11 

l8o 

ll*  .56 

0, 

0 

17.7 

28 

180 

ll*  .56 

1. 

00 

19.0 

U2 

90 

15. W 

0. 

0 

19.3 

12 

180 

15.1*8 

0. 

0 

19.1 

29 

180 

15.1*8 

1. 

00 

21.5 

1*3 

90 

16.1*0 

0. 

0 

21.3 

13 

l8o 

16.1*0 

0. 

0 

22.3 

30 

180 

16.1*0 

1. 

00 

23.2 

1*1* 

90 

17.31 

0. 

0 

21*.? 

11* 

180 

17.31 

2l*.0 

1*5 

90 

18.18 

?6.o 

15 

180 

18. ’8 

25.8 

1*6 

90 

19.00 

28.3 

16 

180 

19.06 

29.0 

1*7 

90 

19.93 

30.3 

17 

180 

19.93 

31.0 

1.8 

90 

20.78 

32.3 

19 

180 

0.78 

32.2 

1*9 

90 

21.61* 

31*.  6 

19 

l8o 

21.61* 

3>*. 3 

50 

90 

22.1.8 

36.3 

20 

l8o 

22.1*8 

36.2 

52 

270 

22.1*8 

36.3 

BB005 

ll* 

1' 

8 

8.8 

6.8 

15 

9.22 

7.1* 

l6 

10.1*5 

9.5 

17 

11.61 

11.7 

18 

12.65 

13.6 

19 

13.65 

15.1* 

20 

ll*.  56 

18.0 

21 

15.!*8 

19.2 

22 

16.1*0 

22.3 

23 

17.31 

23.2 

2l* 

18.18 

26.1* 

25 

19.06 

27.3 

26 

19.93 

30.  ^ 

27 

20.78 

31.1* 

28 

21.61* 

35.1 

29 

22.1*8 

36.7 

• Azimuth  relative  to  point  ir  center  of  large  concrete  block  which  was  base  of  wire-breaker 
displacement-time  stand  used  in  test  shots  BdOOI  and  BB002.  Origin  of  coordinate  system  at 
ground  zero. 

**  The  reference  plane  is  the  horizontal  plane  which  includes  the  point  source  of  the  explosion 
(72  cm  above  original  ground  surface).  Positive  values  indicate  height  above  the  plane.  This 
plane  was  the  same  for  each  test. 


Table  7 

Data  from  Overpressure  (OP)  Gages 


Test 

Shot 

OP  Gage 
Designation 

Azimuth 

deg* 

Distance  from 
Ground  Zero 
m 

Height  of  Gage 
Above  or  Below 
Reference  Plane** 
m 

Activation 
Time  of 
Gage 
msec 

Peak 

Over- 

pressure 

psi 

BB001 

1* 

23 

8 

.8 

0.76 

7.0 

HI 

Of 

1 

90 

8.8 

-0.7>* 

5.7 

2 

180 

8.8 

-0.71* 

6.0 

3 

270 

8.8 

-0.71* 

6.0 

BB002 

1 

93 

5.1*9 

-0.91* 

2.5 

2 

183 

5.1*9 

-0.91* 

2.5 

3 

23 

8 

.8 

0.56 

6.6 

BBOO  3 

1 

93 

5.1*9 

-0.60 

2.9 

I! 

29 

o 

183 

-0.60 

2.1* 

139 

3 

273 

0 

.0 

2.8 

121* 

BBOOU 

1 

138 

0 

.0 

2.7 

128 

2 

183 

0 

.0 

2.8 

ll*2 

3 

273 

0 

.0 

2.8 

126 

BB005 

1 

108 

-0.60 

2.7 

122 

2 

293 

-0.60 

2.7 

127 

3 

207 

8.8 

-0.60 

6.5 

1*6 

5 

331 

8, 

.8 

0, 

.0 

6.1* 

37 

It 

3l*l* 

8 

.8 

6.1* 

31* 

6 

303 

12 

.2 

12.1 

16 

BB006 

1 

108 

5.U9 

3.0 

133 

2 

292 

5.1*9 

2.9 

11*0 

3 

207 

8, 

.8 

6.1* 

52 

BB007 

1 

108 

5.65 

2.7 

ND 

3 

207 

10, 

.7 

9.7 

32 

2 

292 

6, 

.19 

3.3 

100 

* Azimuth  relative  to  point  in  center  of  large  concrete  block  which  was  base  of 
wire-breaker  displacement-time  stand  U3ed  in  test  shots  BB001  and  BBOO?.  Origin 
of  coordinate  system  at  ground  sc-ro. 

**  The  reference  plane  is  the  horizontal  plane  which  includes  the  point  source  of 
the  explosion  (72  cm  above  original  ground  surface).  Positive  values  indicate 
height  above  the  plane  and  negative  below  the  plane.  This  plane  was  the  same  for 
each  test. 

+ ’!o  data. 


Table  8 

Time-of -Arrival  (TOA)  Data  Used  to 
Evaluate  Shock-Front  Symmetry 


Test 

Shot 

TOA  Gage 
Designation 

Azimuth 

deg* 

Distance  From 
Ground  Zero 

m 

Arrival 

Time 

msec 

BB001 

9 

0 

8, 

,8 

6.1 

10 

90 

6.2 

1 

180 

6.2 

8 

270 

6.0 

BB002 

53 

0 

6.7 

51 

270 

6.7 

1*1 

93 

ll*. 

.56 

17.2 

11 

183 

ll*. 

.56 

17.7 

U2 

90 

15. 

,U8 

19.3 

12 

180 

15. 

CO 

-X 

19.1 

U3 

90 

16. 

.1*0 

21.3 

13 

180 

16. 

,1*0 

22.3 

1*1* 

90 

17. 

,31 

2U.2 

lU 

180 

17. 

,31 

21*.  0 

U5 

90 

18. 

,18 

26.0 

15 

l8o 

18. 

,18 

25.8 

1*6 

90 

19. 

,06 

28.3 

16 

180 

19. 

06 

29.0 

1*7 

90 

19. 

93 

30.3 

17 

180 

19. 

93 

31.0 

1*8 

90 

20. 

78 

32.3 

18 

l8o 

20. 

78 

32.2 

1*9 

90 

21. 

61* 

34.6 

19 

180 

21. 

61* 

34.3 

50 

90 

22. 

1*8 

36.3 

20 

180 

22. 

CO 

36.2 

* Azimuth  relative  to  point  in  center  of  large  concrete  block  which  was 
base  of  wire-breaker  displacement-time  stand  used  in  test  shots  BB001 
and  BB002.  Origin  of  coordinate  system  at  ground  zero. 


Table  9 

Time-of -Arrival  Data  Used  to  Examine  Inclination  of  Shock  Front 


TOA  Gage  Designation 

Shot 

Lover 

Upper 

BB001 

2 

5 

3 

6 

U 

7 

BB002 

11 

28 

12 

29 

13 

30 

Activation  Time  of  Gage,  msec 


Lover 

Upper 

8.6 

10.2 

12.1 

lU.O 

15.0 

16.8 

17.7 

19.0 

19.1 

21.5 

22.3 

23.2 

Table  10 

Coaparlson  of  WES  and  BRL  Tlac-of-Arrlval  Data 


Test 

Distance  frost 

Shock-Front 

Tiae  of  Arrival,  nsec 

Shot 

Ground  Zero.  ■ 

Measured  by  WES* 

Predicted  from  BRL  Data 

BB001 

8.8 

6.2 

7.3 

8.8 

6.0 

7.3 

8.8 

6.1 

7.3 

8.0 

6.2 

7.3 

10.8 

8.6 

10.6 

11.05 

10.2 

10.9 

12.3 

12.1 

13.5 

12.52 

14.0 

13.8 

13.8 

15.0 

16.5 

14.00 

16.8 

16.7 

BB002 

8.8 

6.7 

7.3 

8.8 

6.7 

7.3 

10.1*9 

11.4 

10.3 

11.65 

13.1 

12.3 

12.68 

15.2 

14.2 

13.68 

17.4 

16.2 

11*.  56 

17.2 

18.0 

14.56 

17.7 

18.0 

14.57 

19.0 

18.1 

15.1*8 

19.1 

20.0 

15.1*8 

19.3 

20.0 

15.51 

21.5 

20.0 

16.1*0 

21.3 

22.4 

16.1*0 

22.3 

22.4 

16.1*3 

23.2 

22.4 

17.31 

24.0 

24.4 

17.31 

24.2 

24.4 

18.18 

25.8 

26.5 

18.18 

26.0 

26.5 

19.06 

28.3 

29.0 

19.06 

29.0 

29.0 

19.93 

30.3 

30.0 

19.93 

31.0 

30.0 

20.78 

32.2 

32.8 

20.78 

32.3 

32.8 

21.61* 

34.3 

35.0 

21.61* 

34.6 

35.0 

22.1*8 

36.2 

37.5 

22.1*8 

36.3 

37.5 

22.1*8 

36.3 

37.5 

BB005 

8.80 

6.8 

7.3 

9.22 

7.4 

8.0 

10.1*5 

9.5 

10.0 

11. '.1 

11.7 

12.0 

12.65 

13.6 

14.0 

13.65 

15.4 

16.0 

14.65 

18.0 

18.0 

15.48 

19.2 

20.0 

16. 40 

22.3 

22.4 

17.31 

23.2 

24.4 

18.18 

26.4 

26.5 

19.06 

27.3 

29.0 

19.93 

30.2 

30.0 

20.78 

31.4 

32.8 

21.64 

35.1 

35.0 

22.48 

36.7 

37.5 

• Fr  «n 

Table  6. 

Table  11 

WES  and  BRL  Peak  Overpressure  Data 


Test  Distance  from  Peak  Overpressure,  psi 


Shot 

Ground  Zero,  m 

Measured  by  WES* 

Predicted  from  BRL  Data 

BB003 

5.1*9 

12l* 

125 

5.52 

129 

123 

5.1*9 

139 

125 

BBOOl* 

5. »*9 

126 

125 

5.1*9 

128 

125 

5.1*9 

ll*2 

125 

BB005 

5.52 

122 

123 

5.52 

127 

123 

8.80 

31* 

1*3 

8.80 

37 

1*3 

8.82 

1*6 

1*3 

12.20 

16 

21 

BB006 

5.1*9 

133 

125 

5.U9 

ll»0 

125 

8.80 

52 

1*3 

BB007 

6.19 

100 

97 

10.7 

32 

28 

• From  Table  7. 
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* Log  mass  including  mounted  instrumentation. 

••  ACC  * acceleration  gage;  VEL  = velocity  gage;  DSPMT  = displacement-time 
t Smooth  log  (turned  down  on  lathe  to  remove  log  hark). 


Table  13 


Comparison  of  Horizontal  Velocities  from  Left  and  Right  Gage 
Positions  Measured  vith  Displacement-Time  Gages 


Log  Mass 

Distance  of  Log 
from  GZ 

Final  Velocity 
n/sec 

..  g 

■ 

Left 

Right 

9,535 

k.k 

23.9 

25.7 

9,520 

8.8 

7.7 

6.9 

28,367 

8.8 

U.l 

U.l 

32,590 

8.8 

3.7 

3.1< 

38,206 

12.2 

1.8 

2.3 

APPENDIX  A:  TREE  STEM  FAILURE 


1.  This  appendix  describes  the  development  of  a failure  criterion 
for  a tree  stem  under  dynamic  loading  by  a shock  front.  The  failure 
criterion  is  based  on  the  mechanical  properties  of  the  tree  stem  and 
does  not  specify  the  height  on  the  stem  at  which  failure  occurs.  Esti- 
mation of  this  height  requires  the  solution  of  a fourth-order  partial 
differential  equation  (PDE)  that  describes  the  deformation  of  the  stem 
under  loading;  the  solution  of  this  PDE  using  the  finite  difference 
method  is  demonstrated.  The  solution  of  the  PDE  requires  a theoretical 
loading  function  that  simulates  the  effect  of  the  shock  front  on  the 
stem;  the  loading  function  is  derived. 

Failure  Criterion 


2.  The  failure  curvature  of  a tree  stem  under  dynamic  loading 
can  be  derived  by  treating  the  tree  as  a right-circular  cylindrical 
wood  beam.  The  bending  moment  at  any  cross  section  of  a beam  is  ex- 
pressed in  equation  form  as 


“b 


O x I 

r 


(Al) 


where 

* static  bending  moment 

0 = maximum  stress  occurring  at  a section 

1 * irr  /U  * moment  of  inertia  of  circular  cross  section  of  a 

beam  with  respect  to  the  neutral  axis 

r = radius  of  beam  (distance  from  the  neutral  axis  to  the  outer- 
most fiber) 


Substituting  the  modulus  of  rupture,  S , for  the  maximum  stress  in 
Equation  Al  yields 

..  S x I 

c r (A2) 


Al 


where  Mc  « bending  moment  required  for  failure.  A mathematical  descrip- 
tion of  stem  curvature  in  response  to  impulse  loading  is  a prerequisite 
to  the  analysis  of  stem  failure.  The  bending  moment  at  any  height  above 
the  ground  for  any  time  is  described  by 


(A3) 


where 


E = Young's  modulus  parallel  to  longitudinal  axis  of  beam 
2 2 

o y/3x  * geometric  curvature  of  beam 

y(x,t)  * deflection  perpendicular  to  the  longitudinal  axis  of  the 
beam 

x * distance  along  longitudinal  axis  of  beam  from  base  of 
beam  (origin  of  coordinate  system) 

t * time  after  initial  application  of  driving  force  (shock 
front ) 

By  introducing  Equation  A2  into  Equation  A3  and  rearranging  terms,  the 
equation  for  beam  (stem)  failure  curvature  Kc  is 


K 


c 


S_ 

Er 


(AU) 


2 2 

Whenever  the  value  for  the  term  3 y/3x  > Kc  , stem  failure  will  occur 
for  given  values  of  S , E , and  r . Equation  A1*  is  time  independent 
and  does  not  specify  a unique  height  on  the  stem  at  which  failure  occurs; 
also  Equation  AU  is  not  a function  of  the  dynamic  loading.  These  re- 
strictions should  be  removed  so  that  Equation  AU  can  be  used  to  predict 
stem-remnant  height  after  detonation  of  a munition. 

Partial  Differential  Equation 


3.  The  geometric  curvature  of  a beam  can  be  explicitly  related 
to  the  dynamic  loading  of  the  beam  for  unique  values  of  x and  t by 
consideration  of  the  dynamic  equilibrium  equation  of  the  beam: 

+ iL  F(x.t)  (A5) 


A2 


where 


o * pA  * mass  per  unit  length  of  cylinder 
p * density  of  the  wood  beam 
A = cross-sectional  area  of  the  beam 

3y(x,t)/3t  * velocity  of  stem  perpendicular  to  the  longitudinal 
axis  of  the  beam 

F(x,t)  = function  that  characterizes  the  dynamic  loading  of 
the  stem  by  the  shock  front.  F(x,t)  has  units  of 
force  per  unit  length 

The  first  term  of  Equation  A5  represents  the  inertial  force  per  unit 
length  required  to  accelerate  an  elemental  mass  (i.e.  the  mass  associ- 
ated with  a unit  length)  in  a horizontal  direction  (y)  for  a given  unit 
time.  The  second  term  is  the  second  derivative  of  the  bending  moment 
(Equation  A3).  This  terra  represents  the  internal  load  intensity  for  a 
given  beam  cross  section.  In  the  static  condition,  the  vector  resul- 
tant of  all  loads  is  zero,  i.e.  IF  * 0 . The  system  is  in  equilibrium 
because  there  are  no  unbalanced  forces.  Any  excess  in  the  external 
load  over  the  internal  load  must  be  balanced  by  force  that  produces 
acceleration  and  bending  of  the  beam.  By  taking  the  indicated  partial 
derivatives  in  Equation  A5, 


If  p and  E are  assumed  to  be  constant  and  the  cross-sectional  area 
2 

of  the  stem  is  nr  , where  r is  a function  of  x only,  then  the 
following  working  form  of  the  PDE  is  generated: 


TOTAL.  EFFECTIVE  PRESSURE  P. 


By  representing  the  tree  stem  as  a right-circular  cylindrical  bean,  the 
moment  of  inertia  is  constant  for  all  x . Therefore,  the  second  and 
third  terns  of  Equation  A6  go  to  zero,  and  the  equation  becomes 


2 

p«r 


|3v(x,t)"|  _ 

aSrU.t) 

L at  J * EI 

. 1* 
L 3x 

F(x,t) 


(A7) 


Loading  Function 


U.  The  solution  of  Equation  A7  requires  a mathematical  expres- 
sion for  F(x,t)  that  describes  the  loading  of  the  beam  by  the  shock 
front.  F(x,t)  , hereafter  referred  to  as  the  loading  function,  can  be 
established  in  terms  of  the  total  effective  pressure  impulse  (JT)  as 
follows : 

T 

,P 

\ JT  * / P(t)  dt 

o 


where 

positive  pulse  duration 

total  effective  pressure  incident  upon  the  beam 
By  assuming  that  P(t)  rises  to  a maximum  value  instantaneously,  the 
impulse  is  the  area  of  simple  triangle  (Figure  Al).  Therefore, 


T 

P 

P(t) 


Figure  Al.  Triangular  assumption 
used  to  calculate  the  total  ef- 
fective pressure  impulse  (J^,) 


Ah 


TIME  I ♦! 


PT  \Tp 
max  r 


Solving  for 

max 


max 


T 

P 


Multiplying  this  equation  by  the  diameter  of  the  stem  (2r)  produces  an 
expression  with  terms  of  force  per  unit  length,  which  are  the  correct 
units  for  Equation  A7*  Then 

F(x,t)T  = (2r)/P  \ = (A8) 

max  \ max/  p 

The  value  computed  for  Equation  A8  is  the  maximum  force  per  unit  length 

experienced  by  the  stem  during  the  positive  pulse  duration  (T^). 

F(x,t)T  is  a constant  at  a given  point  on  the  stem  for  the  positive 
max 

pulse  duration.  It  is,  therefore,  necessary  to  construct  a loading  func- 
tion that  includes  F(x,t)T  , but  is  a function  of  time. 

max 


5.  The  loading  function  F(x,t)  was  constructed  in  two  com- 
ponents. The  first  component  describes  the  effect  of  the  shock  front 
on  the  stem  from  the  time  the  shock  front  first  comes  in  contact  with 
the  stem  until  the  amplitude  of  the  shock  front  pulse  reaches  a maxi- 
mum value  F(x,t),j,  . This  time  interval  is  defined  as  the  rise 

max 

time  (Tr).  The  second  component  describes  the  effect  of  the  shock 
front  on  the  stem  from  the  time  when  the  amplitude  of  the  shock- 
front  pulse  is  a maximum  until  the  time  when  the  pulse  first  decays 
to  zero.  This  time  interval  is  defined  as  the  decay  time.  The 
time  period  from  initial  contact  of  the  stem  by  the  shock  front 
until  the  shock  front  passes  its  peak  amplitude  and  first  returns 
to  zero  is  defined  as  the  positive  pulse  duration  (T^).  The  above 
requirements  can  be  expressed  mathematically  as: 


A5 


At  t = 0 


At  t * T 

r 

At  t ■ T 

P 


F(x,t)  = 0 
F(x,t)  = F(x,t)T 

max 

F(x,t)  = 0 


(A9) 


6.  It  is  necessary  to  derive  appropriate  scale  factors  which, 

when  multiplied  by  F(x,t)T  , will  satisfy  the  requirements  of 

max 

Equations  A9.  The  scale  factors  must  characterize  the  physical  phenom- 
ena of  the  rise  and  decay  components  of  the  loading  function.  The  load- 
ing function  can  then  be  expressed  mathematically  as: 


F(x,t) 


rise 


F(x,t)T  SFr 
max 


F(x,t) 


decay 


F(x,t)T  SFd 
max 


where 

SFn  = rise  scale  factor 

n 

SFp  * decay  scale  factor 

7.  To  satisfy  Equations  A9,  the  following  restrictions  must  be 
placed  on  the  scale  factors: 


At  t = 0 
At  t = T 

r 

At  t = T 

P 


sfr  = o 
sfr  = sfd  * 


(A10) 


The  value  of  SF_  was  hypothesized  to  increase  as  the  stem  loading  area 
n 

increased.  SF^  was  assumed  to  decrease  linearly. 

8.  The  stem  loading  area  is  assumed  to  be  a plane  perpendicular 
to  the  horizontal  component  of  the  total  effective  pressure  impulse  and 
confined  to  the  chord  of  the  stem  surface  and  length  differential  (dx) 
of  the  stem  (Figure  A2) . The  rise  time  Tp  is  assumed  to  be  the  ratio 
of  the  stem  radius  to  the  shock- front  velocity.  The  time  required  to 


A6 


Figure  A2.  Calculation  of  time  required  to  load  a 
given  differential  stem  area 

load  a given  area,  , is  calculated  as  follows: 

T = It2  - (T  - t)2 
a -W  r r 

where 

= time  required  to  load  differential  stem  area 
Tr  = rise  time  of  positive  pulse 
Simplifying, 


The  dimensionless  square  root  is  equivalent  to  SFD  because  it  obeys 

r\ 

the  required  conditions  of  Equation  A10  and  increases  as  the  loading 
area  increases.  Therefore, 


F(x,t) 


rise 


max 


(All) 


9.  The  decay  component  of  the  loading  function  can  be  derived 

from  the  equation  of  the  line  from  fr,  ,F(x,t)^,  "jtoCT  ,0)  (see 

L maxJ  F 

Figure  A3).  The  slope  of  the  line  is 


A7 


Figure  A3.  Graphic  representation  of  loading 
function 


m = - 


max 


T - T 
P r 


The  intercept  is 


B - Ft  + l - Ft  + m x Tr 
max  max 


Therefore, 
F(x,t) 


max 


decay  mt  + B " " \Tp  - Tr/ 


t + 


p i — E22 — i t 

T VT  - T / r 
max\  p r/ 


max 


t - T \ 

rrf  for  T < t < T 
p r/  y 


10.  Equations  All  and  A12  satisfy  the  requirements  in  para- 
graphs 5 and  7,  above,  but  the  derivatives  of  the  functions  are  di; 

continuous  at  t = T , i.e. 

r ’ 


3F 


rise 


at 


decay 

3t 


A8 


( A12 ) 


(A13) 


Physically,  this  means  that  at  t * T^  the  rise  and  decay  functions 
are  not  changing  at  the  same  rate  with  respect  to  time.  This  problem 
was  rectified  by  multiplying  the  rise  component  of  the  loading  function 
by  the  decay  scale  factor  (SF^).  The  rise  component  of  the  loading 
function  then  becomes 


F(x,t) 


rise 


for  0 < t < T (AH*) 

r 


Equation  Al1*  then  satisfies  the  equality  (Equation  A13).  SF^  is  very 
close  to  unity  between  t = 0 and  t = Tf  ; therefore,  it  is  numerically 
insignificant  in  the  solution,  but  it  does  cause  the  peak  value  of 
loading  to  occur  slightly  before  t = Tr  (Figure  A3). 

11.  Two  further  refinements  were  made.  First,  the  total  effec- 
tive pressure  impulse  discussed  thus  far  has  assumed  radial  symmetry 
with  respect  to  the  source  of  the  shock  front  (center  of  explosion). 

At  heights  on  the  stem  different  from  the  height  of  the  center  of  the 
explosion,  error  will  be  introduced  into  the  computations  if  the  radial 
value  of  J,p  is  used;  therefore,  to  eliminate  this  error,  the  horizon- 
tal component  of  must  be  substituted  into  the  loading  function,  i.e. 


JT  = JT  COS  4 

horizontal 


where  4 is  the  angle  between  the  horizontal  and  the  point  of  interest 
on  the  tree.  Second,  the  value  of  J^,  is  not  constant  for  all  points 
on  the  stem  because  the  shock  front  decays  as  the  reciprocal  of  tL° 
square  of  the  propagation  distance;  and  since  the  distances  from  the 
center  of  the  explosion  to  any  two  points  on  the  stem  are  different, 

JT  must  also  be  different.  Therefore, 


TL 


A9 


where 


JTL  * effective  total  pressure  impulse  at  any  point  on  the  stem 

J,pp  * effective  total  pressure  impulse  at  a height  on  the  stem 
equal  to  the  height  of  the  center  of  the  explosion  (EC) 

R = horizontal  distance  from  EC  to  tree  stem 

L = radial  distance  from  EC  to  point  of  interest  on  tree  stem 


These  two  modifications  are  shown  in  detail  on  pages  16  and  17  of 
Reference  1.* 

12.  'Hie  resulting  loading  function  is 


F(x,t)  . 

’ rise 


^rJ^R  cos  <p 


T L 
P 


F(x,t) 


decay 


UrJTRR  cos  4* 

T L2 
P 


(-£*) 


for  0 < t < T 


for  T < t < T 
r p 


A value  for  JTR  was  calculated  by  the  following  equation: 

JTR  ■ JCd  + Pr  (r) 

where 

J = dynamic  pressure  impulse 

C,  = drag  coefficient  of  stem 
d 

P = reflected  pressure 
r 

U = shock-front  velocity 

All  values  apply  to  the  point  on  the  stem  that  equals  the  height  of  the 
EC.  JC^  results  from  the  air  drag  on  the  stem  within  the  shock  front. 
The  second  term  in  the  equation  is  due  to  the  pressure  differential  be- 
tween the  side  facing  the  incident  shock  front  and  the  lee  side  of  the 
tree  stem.  The  term  2r/U  is  the  shock-front  transit  time  across  the 
stem  (see  Reference  1,  pages  25  and  26). 


* References  in  this  appendix  are  cited  at  the  end  of  the  main  text. 
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Finite  Difference  Solution 


13.  The  fourth-order  differential  equation  describing  the  dynamic 
deformation  of  a right-circular  cylindrical  beam  (Equation  A7 ) is  dif- 
ficult to  solve  in  closed  form  for  an  arbitrary  function  F(x,t)  . 
Therefore,  a numerical  scheme  was  employed  to  obtain  a solution  of 
Equation  A7.  The  derivatives  were  represented  in  finite  difference 
form  as  follows: 


pirr 


G“ 


- 2Y  + Y 

-i,j  fu  Vij 


At 


) 


* * *1.1  - ‘"l.J*!  * W.  . p (A15) 

^ V ^ 5) 

where 

Y^  = horizontal  displacement  at  time  step  i and  nodal  point 
j above  the  ground 

i = time  subscript 

j = distance  subscript 

At  = time  increment 

AH  - distance  increment  along  the  beam 

F.  = loading  function  in  terms  of  discrete  values  (F(x,t) 
implies  a continuous  function) 

Note  that  upper  case  Y's  are  used  for  the  finite  difference  solution  and 
lower  case  y's  for  differential  representation. 


Boundary  Conditions 


ll*.  Associated  with  the  solution  of  the  fourth-order  PDE  are  four 
boundary  conditions  related  to  the  independent  variable  x and  two  re- 
lated to  the  independent  variable  t (more  properly  an  initial  value 
problem).  Consider  the  stem  model  as  in  Figure  kk.  The  stem  is  free 
to  rotate  at  the  base,  but  not  to  translate.  The  following  boundary 
conditions  can  be  deduced  from  the  diagram:  (a)  At  x = 0 , y = 0 , 
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Figure  kh.  Stem  model 


(b)  the  bending  moment  at  x = 0 (base  of  stem)  is  zero,  (c)  the  bend- 
ing moment  at  x = L (top  of  stem)  is  zero,  and  (d)  the  shear  at  x = L 
is  zero.  These  can  be  expressed  mathematically  as 


At  x = 0 , y(0,t)  = 0 and  El(x)  = 0 

3x 


At  x = L , El(x)  = 0 

3xJ 


(Al6a) 


and 


ei(x)  a3y^».^  = 0 

The  statement  (Equation  Al6b)  is  true  because: 


(Al6b) 
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2 2 

But  at  x = L , 3 y/3x  = 0 . Therefore,  expression  Al6b  is  valid. 
Applying  finite  difference  notation  to  the  boundary  conditions: 

At  x = 0 (base  of  stem,  i.e.,  J = 2) 


and 


Y - 2Y  + Y 

-l.li.-J: _lli 

AH2 


0 


Therefore, 

Y = -Y 

i,J-l  xi,J+l 


At  x * L (top  of  stem,  i.e.,  J = n , where  n is  the  last  node 
point) 


r - 2Y  + Y 

Y1 J--  Y-U±l.  . 0 a:d 


AH 


-Y  + 2Y  - 2Y  + Y 

i,J-2  U3+1 Ilii+2  . 0 

2AHJ 


(A17) 


Then, 


Y.  = 2Y.  . - Y.  . , 
i,j+l  i,j  i,J-l 


and  from  Equation  A17 


Y.  4.0  = Uy.  , - UY.  . - + Y.  . . 
i , J+2  i»j  i»J-l 


15.  The  initial  value  problem  relevant  to  time  is  solved  by 
considering  that  the  displacement  and  velocity  at  t = 0 are  zero. 
These  can  be  expressed  mathematically  as 

Y(x,0)  = 0 and  Y(x,t)  = 0 for  t < 0 


and 

3Y(x.t ) _ 
3t 
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Hence,  at  t = 0 , i.e.,  i=l,  the  initial  conditions  can  be  ex- 
pressed in  finite  difference  form  by 


ij 
i-U  = 


and 

-Y  + Y 

■>i.J = o 

2At 


Therefore, 


Y = Y 

i+l,J  i-l.J 


Sample  Solution  of  Finite  Difference  Equation 


16.  The  six  boundary  conditions,  plus  the  finite  difference  form 
for  the  PDE,  are  sufficient  to  generate  values  for  the  deflection  and 
its  derivatives  for  all  times  greater  than  t = 0 at  selected  heights 
up  the  tree.  The  sequence  of  mathematical  operations  is  to  calculate 
the  deflection  at  each  of  the  node  points  for  a given  time,  i.e.,  the 
j subscript  is  incremented  from  J = 1 to  j = n for  a given  value 
of  i . Solving  for  the  acceleration  term  from  Equation  A15: 


i.J 


* i.J 
2 

pirr 


Er 

P 


P 


- i+Y 


luhL 


6Yi  - + Yj . j+2 


AH 


‘) 


Also  from  finite  difference  notation: 


Y - 2Y  + Y 

J 

At2 


Therefore, 


ri«.j = + 2Iu 


- Y 


i-l.J 
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Thus,  the  deflection  in  the  next  time  increment  can  be  generated.  Con- 
sider the  following  graphic  representation  (Figure  A5)  of  eight  node 
points  prior  to  the  initial  calculation.  The  initial  deflections 
(at  t * 0)  are  assumed  to  be  zero,  i.e.,  Y.  . , where  J varies  from 
3 to  10.  From  the  boundary  conditions  Y^  ^ * ^1  11  * an<*  l are 

also  zero.  The  deflections  before  time  zero  (Y  .)  are  also  assumed  to 

o,  j 

be  zero.  The  time  zero  here  specifically  refers  to  the  time  at  which 
the  shock  wave  arrives  at  a node  point  on  the  stem  and  differs  from 
the  detonation  time  for  the  bomb  by  the  amount  of  time  required  for 
the  shock  front  to  travel  from  ground  zero  (GZ)  to  the  node  point  in 
question. 


17. 


Since  all  the  Y,  ,'s 
* 9 J 


are  zero  for  the  initial  calculations. 


Figure  A5.  Graphic  representation  of  initial 
data  at  t = 0 . As  discussed  in  the  text, 
all  Y.  . values  are  zero 

i .J 
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F 

Y .-ill 

i,J  2 

pvr 


where  F.  . is  different  for  each  node  point.  The  difference  in  the 
* * J 

driving  function  F.  . at  each  node  point  arises  from  the  fact  that 

i » j 

the  distance  that  the  shock  front  must  travel  is  different  at  each 
point,  and  the  angle  between  the  horizontal  and  incident  front  is  dif- 
ferent for  each  node  point.  The  following  calculation  can  now  be  made 


2 

pirr 


2 2 

where  B = At  /pur  . The  deflections  for  the  first  two  time  steps 

are  shown  in  Figure  A6.  The  data  enclosed  by  ellipses  represent  the 

data  necessary  to  obtain  the  deflection  for  Y_  _ . The  five  points 

*•»  ( 

in  the  vertical  ellipse  \ ield  the  acceleration;  then  using  a backward 
difference,  the  actual  deflection  is  calculated.  The  specific  details 
of  obtaining  Y2  ^ are  as  follows: 


+ 6Y 


P 


AH 


But 


Yl,5  = *1,6  = *1,7  ” *1,8  = *1,9 
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Therefore , 


Then 


*bW 
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pwr 


pirr 


+ 2(0)  - (0)  = BF1>7 
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Figure  A6.  Calculations  for  displacements  at  t * At^  . 
Numbers  enclosed  in  elliptical  areas  show  data  used  to 

calculate  Y„  7 

^9  I 

2 2 

where  3 * At  /pur  . The  four  derivatives  of  the  deflection  can  then 
be  calculated  by  using  central  differences.  The  second,  third,  and 
fourth  derivatives  correspond  to  the  bending  moment,  shear,  and  load, 
respectively,  provided  the  central  difference  is  multiplied  by  the 
flexural  rigidity  (El).  The  calculations  continue  in  this  manner  until 
the  limit  of  the  time  interval  of  interest  is  reached. 

Stability  of  Finite  Difference  Equation 

18.  The  stability  of  the  finite  difference  equation  for  the 
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given  selected  differences  can  be  examined  using  the  von  Neuman  neces- 
sary condition  for  stability  of  thin  elastic  beams  from  Reference  8; 


aAt 


1/2 


where 

a = EI/u 

U = mass  density  per  unit  length 
At  = time  increment 
AH  = distance  increment 

Then 


For  the  reference  case  solution  (see  Reference  l)  used  to  generate 
failure  heights  versus  distance  from  ground  zero  (DGZ), 


E * 1.2  x 10 


11 


p * 0.8 


Then, 


(1.2  x lOU)r2 


< 1/2 


Therefore,  the  stability  condition  for  this  problem  is 


< 0,258  * 10"5 


(A18) 


Output  from  the  computer  solution  of  Equation  A15  was  verified  using 
Equation  Al8.  Violation  of  Equation  Al8  yielded  deflections  that  were 
meaningless . 


Al8 


APPENDIX  B:  NOMOGRAPHIC  SOLUTION  FOR 

PREDICTION  OF  STEM-REMNANT  HEIGHT 


1.  This  appendix  describes  the  development  of  a nomographic 
solution  for  the  prediction  of  tree  stem-remnant  height  for  wide  varia- 
tions in  weapon  yield  and  stem  properties.  This  solution  uses  the  stem 
failure  criterion  and  the  finite  difference  solution  of  the  partial  dif- 
ferential equation  (PDE)  describing  stem  motion  developed  in  Appendix  A. 

2.  The  finite  difference  solution  for  the  PDE  is  a lengthy  itera- 
tive technique  that  is  amenable  to  a computer  solution  in  lieu  of  manual 
calculations.  The  resulting  computer  output  yields  the  deflection  of 
the  stem  and  its  derivatives  at  selected  heights  up  the  tree  for  a given 
time.  The  second  derivative  of  the  deflection  is  the  curvature;  hence, 
data  for  curvature  times  stem  radius  versus  height  on  the  stem  at  which 
the  curvature  occurs  are  available  for  a given  stem  diameter.  Using  the 
output,  a family  of  curves  was  generated  for  selected  stem  diameters 
relating  maximum  curvature  multiplied  by  stem  radius  versus  height  on 
the  stem  at  which  the  maximum  curvature  occurred  (Figure  Bl). 

3.  The  curves  for  Figure  Bl  for  a given  stem  diameter  were  ob- 
tained as  follows:  For  each  time  increment  (At),  the  computer  print- 

out provided  the  curvature  for  each  node  point  on  the  stem.  The  node 
at  which  the  maximum  curvature  occurs  and  the  corresponding  height  of 
the  node  provide  the  data  point  for  that  At. 

*4.  The  model  for  the  stem  consists  of  a right-circular  cylinder 
constrained  at  the  base  by  a single  point,  such  that  the  base  is  free 
to  rotate  but  not  translate  (Appendix  A).  Therefore,  as  the  shock  front 
interacts  with  the  model,  the  point  of  maximum  curvature  must  move  up 
the  stem.  This  phenomenon  is  similar  to  whiplash.  Examination  of 
sequential  time  increments  of  the  computer  printout  indicated  that  the 
point  of  maximum  curvature  did  always  move  up  the  stem;  but  because  of 
the  complex  oscillatory  motion  of  the  stem,  the  value  of  the  maximum 
curvature  did  not  always  increase.  Data  points  that  had  a lower  value 
of  maximum  curvature  than  the  previous  time  step,  but  a higher  node 
point,  were  rejected. 


Bl 


HEIGHT  ON  STEM,  M 


Figure  Bl.  Maximum  stem  curvature  multiplied  by  stem  ra- 
dius versus  height  at  which  maximum  curvature  occurred 

5.  The  height  on  the  stem  at  which  the  maximum  curvature  occurs 
and  the  stem-remnant  height  are  considered  to  be  the  same  if  failure 
occurs.  A remnant  height  obtained  from  Figure  Bl  can  be  considered  to 
be  the  remnant  height  associated  with  the  value  of  Kc  x r required  for 
failure  (Appendix  A),  because  the  curves  in  Figure  Bl  represent  maximum 
values  of  curvature  on  the  stem  at  a given  time;  therefore,  if  a calcu- 
lated value  of  the  critical  curvature  times  the  radius  is  entered  on 
the  vertical  scale  of  Figure  Bl  and  has  a corresponding  value  on  the 
horizontal  axis,  it  then  follows  that  this  is  the  remnant  height  asso- 
ciated with  the  calculation  of  K * r because  the  critical  curvature 

c 

and  the  maximum  curvature  are  the  same.  If  K « r does  not  have  a 

c 
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value  on  the  horizontal  scale,  the  stem  did  not  fail  and  the  energy 
absorbed  by  the  stem  from  the  shock  front  is  dissipated  as  vibratory 
motion  or  structural  damage.  In  view  of  the  above  interpretation,  the 
axes  of  Figure  B1  can  be  relabeled  critical  curvature  times  stem  radius 
versus  remnant  height. 

6.  The  curves  generated  in  Figure  B1  by  inspection  of  the  print- 
outs for  each  radius  (r)  are  valid  only  for  the  value  of  Young's  modu- 
lus (E),  total  effective  pressure  impulse  ( J^,) , and  stem  density  (p) 
used  in  the  finite  difference  solution  of  the  PDE.  The  solution  of  the 
PDE  used  to  proauce  Figure  B1  is  termed  the  reference  case  solution. 
Numerical  values  used  in  this  solution  are  termed  reference  case  values, 
i.e.,  reference  case  stem  radius,  reference  case  Young's  modulus,  refer- 
ence case  total  effective  pressuie  impulse,  reference  case  stem  density. 
The  reference  case  values  used  in  the  solution  were  r = 6,  8,  12,  20, 

40,  60,  and  75  cm,  E * 1.2  * 1011  dynes/cm2,  J ■ 4.99  x 10**  dyne 

2 3 1 

x msec/cm  ,*  and  p = 0.8  g/cm  . The  curves  in  Figure  B1  are  not  ap- 
plicable if  a remnant  height  prediction  is  desired  using  values  other 
than  those  reference  case  values  used  in  the  solution  that  generated  the 
curves  in  the  figure.  In  order  to  make  Figure  B1  applicable  for  wide 
variations  in  E , J^,  , and  p , two  options  were  available:  (a)  gener- 
ate curves  for  all  desired  values  of  E , , and  p , or  (b)  derive 

scale  functions  such  that  values  of  E , J^,  , and  p other  than  those 
used  in  the  reference  case  PDE  solution  can  be  used  to  predict  remnant 
heights  using  Figure  Bl.  Option  (b)  is  the  logical  choice. 

7.  Use  of  Figure  Bl  for  specific  values  of  E , p , and 
(other  than  those  used  in  the  reference  case  PDE  solution)  requires 
modification  of  the  failure  criterion  to  account  for  the  variations  in 
E , JT  , and  p . The  modified  failure  criterion  is  as  follows: 


(Bl) 


* References  in  this  Appendix  are  cited  at  the  end  of  the  main  text. 
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where 


K = reference  case  critical  curvature  required  for  stem  failure 
co 

r = reference  case  stem  radius 
o 

J^o  = reference  case  total  effective  pressure  impulse 

PQ  = reference  case  stem  density 

= drag  coefficient  (0.75  was  used  in  the  reference  case 
solution) 

The  subscript  o indicates  a parameter  used  in  the  reference  case 
solution.  The  other  terms  in  Equation  B1  are: 

S = modulus  of  rupture 
E = Young's  modulus 
p = stem  density 
W = weapon  yield 
r = stem  radius 

all  of  which  refer  to  the  specific  case  solution.  The  values  that  refer 

to  the  actual  stem  whose  remnant  height  is  required  are  termed  specific 

case  values.  Also,  J_  = known  value  of  the  dynamic  pressure  impulse  at 

b 

a given  distance  from  GZ  for  weapon  yield  Wfi  . The  curves  in  Figure  8 
of  Reference  1*  were  used  to  compute  J.,  for  the  reference  case  curves 
in  Figure  Bl.  The  derivation  of  Equation  B1  is  shown  in  Reference  1, 
pages  25-3^.  Also,  more  information  on  the  reference  and  specific  case 
solutions  can  be  found  in  paragraphs  60-80  of  Rei - -ence  1. 

8.  The  remnant  height  for  a given  stem  can  be  obtained  by  the 
following  procedure: 

a.  Kcq  x rQ  is  computed  from  Equation  Bl. 

b.  The  diameter  curves  in  Figure  Bl  are  reference  case 
values  and  differ  from  the  specific  case  values  by 


(B2) 


Therefore,  dQ  must  be  calculated  for  a given  d . 


* References  in  this  appendix  are  cited  at  the  end  of  the  main  text. 


The  remnant  height  can  then  be  obtained  by  entering  the 
computed  value  of  Kco  x on  the  vertical  axis  of 
Figure  B1  and  projecting  a horizontal  line  from  this 
point  until  the  proper  reference  case  diameter  curve  is 
intersected.  A line  is  then  projected  parallel  to  the 
vertical  axis  from  the  intersection  to  the  horizontal 
axis;  the  intersection  with  the  horizontal  axis  is  the 
predicted  remnant  height. 

9.  The  procedure  outlined  in  paragraph  8 is  not  conducive  to 
rapid  calculation  of  stem-remnant  height  if  calculational  aids  are  not 
available,  i.e.  in  a combat  situation.  Therefore,  in  an  effort  to 
reduce  computational  difficulties,  the  entire  process  was  reduced  to  a 
nomographic  solution  (Figure  B2),  hereafter  referred  to  as -the  nomograph. 

10.  Equation  B1  can  be  rewritten  as 


\ 


c. 


The  initial  procedure  is  to  graphically  solve  for  the  term  in  the  second 
set  of  brackets.  This  was  accomplished  by  generating  diameter  curves 
for  a continuous  variation  of  weapon  yields  (graph  3 of  the  nomograph). 
The  resulting  numerical  scale  common  to  graphs  3 and  U is  the  computed 
numerical  value  for 


0.75  + 


for  given  values  of  d and  W. 

11.  The  next  step  is  to  obtain  the  ratio  of 


'SJ_To_P 
Eo 


to 


0.75  + 


and  then  divide  this  value  by  J_  (W/W-)^^  to 

D D 


obtain  Kcq  * rQ  . The  graphic  division  process  is  performed  using 
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LCULATKD  VALUC  FOR  TNI  SAMPLE  RROGLfM  IS  AO  CM  I 


« NOT  MEN  DRAWN  A OR  TnC  ACAAON  V*LO  (•>  Of 
ACAAAKO  AT  AC  AO  IWO  SEVERAL  RAlRS  OA  COOROMATCS 
> MULTIPLYING  IACH  IV  ^ , (THIS  VALUC  IS  I t» 
m ANCAC  I IS  IN  METRIC  TONS:  tlM  ROUNDS  EQUALS 
MIRLt,  TMC  COOAOINATCS  70S.  t ON  TMC  W ■ 14 

» ••  IS  AON  ■ A CURVE 

K RNISEMT  RILL  ThCR  SC  OEnEAATEO  THROUGH  TmC 


lltrjL  AT  THC  ROtMT  NMCAC  A HORIZONTAL  UMC  TnNOUOM  TmC  OCHRE  D 

distance  from  onouno  zero  cnoucs  thc  curve  eon  tnc  arrnoraiate  weapon 

TICLO  ONAR  A VIATICAL  LINC  THAT  CUTE NOS  ACROSS  TNC  DlAOONAL  LINES  ON 
ONAAM  » (THIS  ORE  RATION  CALL  YlCLO  ESS  ON  TMC  HORIZONTAL  AXIS  OR  SNAP*  < 
AON  TNC  SAMRLE  RNOOLCM  I 
ORA  AW  S, 

STEA  I IA  NO  CUNVC  Ift  RNCSCNT  AON  TnC  Of  WACO  STEM  DIAMETER  141. 

AN  AAAWQAAlATE  CUAVf  CAN  M ARE  EARED  DV  ANtOlNO  HORIZONTAL  COORDINATES 
ASSOCIATED  NITM  VARIOUS  RE  A RON  YIELDS.  f»l  AS  ROLLOVt: 


THE  NUMERICAL  VALUES  AON  TmC  HORIZONTAL  AXIS  ARC  ARCSCNTEO  ON  TmC  HORI- 
ZONTAL AXIS  COMMON  TO  ORAAMf  s ANO  « IA  ONLY  ONC  RCA  RON  YIELD  IS  MINS 
UTILIZED  The  ENTIRE  CURVE  MCCO  MOT  M DRAWN  UNLESS  IT  RILL  Cf  USE O LATER 
<R<4  A)  EQUALS  I TZ  AOR  THE  SAMRLE  AROOLCM  I 

STEA  z AT  The  INTERSECTION  OA  A HORIZONTAL  LINE  TmROuOm  TmC  GIVEN 
ME  A NON  YlCLO  ANO  CORRECT  STEM  OlAMCTCA  CuRVt  CONSTRUCT  A VERTICAL 
LINC  SO  THAT  IT  RASSCS  TmROuOh  THE  DlAOONAL  UNCI  ON  ORARM  « 

PRAAM  < 

STEA  | CONSTRUCT  A Ol AGONAL  LINC  IRARALLCL  RlTH  THOSE  RRf SENT)  AROM 
THE  CALCULATED  VALUE  OA  ■ 10  * ON  THE  LEAT  VERTICAL  AXIS  IThC  CAL- 


S' ON  THE  LEAT  VERTICAL  AXIS  I TNC  CAL- 


CULATED VALUC  FOR  P* • ir*  FOR  The  SAMRLE  ARORLCM  IS  0 SM  I 


STIRS  FROM  TMC  MTCRSECTION  OR  ThC  DlAROMAL  UNC  M STCR  I ANO  ThC 
VCRTlCAL  UNC  RROJCCTCO  AROM  RNARH  » CONSTRUCT  A HORIZONTAL  LMC  MUCH 
RU  MTCRSCCT  TMC  RIGHT  VCRTlCAL  AXIS  OA  ORAM  « (This  VALUE  IS  0 f 11  AOR 


ITER  I ORAR  A DIAGONAL  LMC  IRARALLCL  RlTH  ThOIC  RRCSCnTI  AROM  TMC 
MTCRSCCT  ION  ROMT  ON  THC  RIGHT  VCRTlCAL  AXIS  OR  ORA  AM  4 TO  TmC  AOmT  On 
TMC  LCRT  VCRTlCAL  AXIS  OR  ORAM*  • IRH0SC  VALUC  IS  MALLCR  OV  A FACTOR 
PA  If 

ITCRS  RHCRC  TMC  DIAGONAL  LMC  M STCR  I MTCRSCCT!  THE  VCRTlCAL 
LMC  CXTCH0MG  FROM  CRAAN  t.  CONSTRUCT  A HORIZONTAL  LMC  THROUGH  ORARh  « 
ITHC  VALUC  OR  TmC  VCRTlCAL  AXIS  OA  PRAAM  # IS  MMU  AOR  TmC  IMRLC  AROPlCM  I 
QRAAWS; 

STCR  I IR  MOCURVC  EXISTS  On  GRARh  • AOR  TmC  STCM  OlAMCTCN  <4#»  CALCU- 
LATED FROM  PRAAM  I.  (fT MATE  A AIR  ROMTS  M TnC  ARRRORRlATI  MTIRVAL  FOR 
TM(  DIAMETERS  GIVEN  ANO  SWITCH  THf  ARRROMIMATE  CURVt  THROUGH  THOSC  AOmTS 
CHOSEN 

^TIR 1 AROM  TMC  INTERSECTION  OA  TMl  «f,  VALUC  (HORIZONTAL  LMC 
(XTENOED  AROM  ORARh  Si  ANO  TmC  NCACRCnCC  CASE  STCM  DlAMCTfR  CURVE.  ORAR  A 
LINE  WHICH  it  RCRRCNOlCULAR  TO  THf  HORIZONTAL  AXIS  OA  QRAAw  • THE  TREE 
RE  Ml  ANT  HEIGHT  IS  TmC  VALUE  ASSOCIATED  RlTH  THE  ROMT  AT  WHICH  THE  RCRRCN- 
OlCULAR LINE  CROSSES  TMC  HORIZONTAL  AXIS  (THE  RC SMART  MClOHT  IS  141  FOR  Th« 
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Preceding  page  blank 


tree  remnant  height  after  explosion  of  a given  size  charge  for  a 
meter  and  strength  at  a given  distance  from  GZ 


logarithmic  paper.  Figure  B3  shows  the  two  possible  configurations 
for  graphic  division;  in  both  cases  A is  the  numerator  and  B the 
denominator.  The  ratio  C is  obtained  by  the  following  procedure: 


Figure  B3.  Use  of  logarithmic  paper  for  graphic  division. 
In  both  cases,  C * A/B  and  scale  A has  a numerical 
value  of  10  times  scale  C 


a.  A line  is  projected  from  the  point  on  the  horizontal 

scale  that  is  the  numerical  value  of  the  denominator  (B). 
This  projected  line  is  perpendicular  to  the  horizontal 
scale. 


b.  A line  is  projected  from  the  point  on  the  vertical  scale 
that  represents  the  numerical  value  of  A (this  scale 
has  numerical  values  of  10  times  the  other  vertical 
scale).  This  projected  line  forms  an  angle  of  1*5  deg 
with  the  horizontal  axis. 


£.  From  the  intersection  of  the  two  projected  lines,  a 
line  is  projected  horizontally  to  the  other  vertical 
scale.  The  numerical  value  of  the  ratio  A/B  is  the 
value  on  the  vertical  scale  where  the  projected  line 
crosses  the  scale.  . . 


12.  Using  the  method  in  Figure  B3a, 


fSJ_To_p 

Ep„ 


is  divided  by 


(graph  k of  the  nomograph).  The  numerical  r ?sult 


of  the  division  is  represented  by  the  scale  common  to  graphs  U and  5 on 
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the  nomograph.  The  value  of 


need  be  calculated  only  once  for 


Eo 


each  tree  species. 

13.  Using  the  method  in  Figure  B3b 


’ \Epo/ 


at 


0.75  + 


is  divided  by  Jgi 
The  dynamic  pressure  impulse 


(graph  5 of  the  nomograph)  to  obtain  Kcq  * rQ 


■ -JfT  ■ 


is  obtained  from  graph  2 by 


entering  the  distance  from  ground  zero  (DGZ)  on  the  vertical  scale  and 

using  the  appropriate  yield  curve.  The  scale  common  to  graphs  2 and  5 

is  the  numerical  value  of  the  dynamic  pressure  impulse  (DPI). 

lU . Graph  6 is  the  same  as  Figure  B1  except  the  vertical  scale  is 

now  labeled  K r and  the  horizontal  scale  is  labeled  tree-rennnnt 
co  o 

height  (in  metres).  Thus  the  value  obtained  in  graph  5 can  be  entered 
into  graph  6,  provided  the  proper  reference  stem  diameter  is  known.  The 
reference  diameter  is  obtained  from  Equation  B2.  Graph  1 is  a solution 
of  this  equation  using  the  method  shown  in  Figure  B3b. 

15.  Additional  curves  for  the  graphs  shown  on  the  nomograph  can 
be  generated  as  follows: 

a_.  Diagonals  for  graph  1 can  be  added  for  other  specific 
case  diameters  as  required. 

b.  Curves  for  yields  other  than  those  shown  in  graph  2 can 
be  generated  from  any  one  of  the  curves  shown  by  using 
the  cube-root  scaling  law  for  both  the  DPI  and  the 
DGZ  (Reference  f). 

c_.  Curves  for  graph  3 for  diameters  other  than  those  shown 


can  be  generated  by  solving 


yields  of  0.1  to  100  metric  tons  (Wg  = 3.6  metric  tons) 
for  the  diameter  required. 


BIO 


d.  Diagonals  for  graph  U can  be  added  for  other  specific 
case  tree  species  as  required. 

£.  Diagonals  for  graph  5 can  be  added  as  required  for  the 
graphical  solution  resulting  from  graphs  3 and  k. 

£.  Curves  for  graph  6 can  be  obtained  for  stem  diameters 
other  than  those  shown  by  substituting  the  required 
diameter  in  the  finite  difference  solution  for  the  par- 
tial differential  equation  describing  the  motion  of  a 
st'm  being  acted  upon  by  an  incident  shock  front  (Appen- 
dix A).  The  data  points  for  the  curve  can  then  be  gener- 
ated from  the  computer  output  of  the  finite  difference 
solution  as  previously  described  (paragraph  3). 

16.  All  of  the  curves  can  be  generated  with  a minimum  of  computa- 
tional difficulty,  except  those  in  graph  6.  Some  interpolation  between 
the  diameter  lines  in  graph  6 is  acceptable  as  a last  resort,  but  with 
an  accompanying  loss  in  the  credibility  in  the  results. 

17-  The  complexity  of  the  nomograph  can  be  reduced  for  situations 
in  which  the  yield  or  stem  diameter  is  constant  or  when  only  one  tree 
species  is  to  be  considered.  Consider  the  case  of  a constant  yield. 
Graph  2 can  be  eliminated  by  converting  the  scale  common  to  graphs  2 and 
5 from  DPI  to  DGZ.  This  can  be  accomplished  by  replotting  values  of  DGZ 
for  a given  yield  curve  along  the  horizontal  axis.  For  constant  yield, 
graph  3 can  also  be  eliminated  by  converting  the  scale  common  to  graphs 


3 and  U from 


to  diameter  values.  This  can  be 


accomplished  by  replotting  the  horizontal  scale  with  selected  diameter 
values.  Similar  schemes  can  be  used  for  situations  in  which  only  one 
stem  diameter  or  species  is  of  interest. 

18.  The  remnant  heights  generated  from  the  nomograph  should  be 
considered  valid  data  only  if: 

a.  The  center  of  the  explosion  is  well  below  the  bottom  of 
the  forest  canopy  and  is  less  than  5 m above  the  ground. 

b.  The  overpressure  is  never  negrtive. 

£.  The  terrain  around  ground  zero  is  flat. 


Bll 


These  constraints  are  explained  in  Reference  1.  In  addition,  the  nomo- 
graph will  not  predict  soil  or  root  failures. 
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Zb  accordant  with  Bt  70-2-3.  par^nqph  6c(l)(b), 
dated  15  Ftoity  1973.  a facsimile  catalog  card 
in  Library  of  Congress  format  la  rnrodneed  below. 
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